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1.  INTRODUCTION 


An  increasing  interest  in  high  speed  towlines  and  towed  devices  by  the 
naval  establishment  has  led  to  the  development  of  low  drag,  faired  (stream- 
lined) towing  systems.  It  is  clearly  important  that  these,  like  all  other 
towing  cables,  must  maintain  a specific  and  predictable  orientation  through- 
out their  operational  speed  ranges. 

Unfortunately,  some  faired  cables  that  have  been  put  into  service  have 
exhibited  dynamic  instabilities.  In  particular,  they  have  experienced  a 
lateral  instability  which  is  commonly  referred  to  as  "kiting".  Although 
steps  to  correct  kiting  have  been  taken  in  the  field, based  on  "engineering 
intuition"  and  experience,  it  is  more  desirable  to  do  so  using  knowledge 
gained  from  studying  the  interactions  between  hydrodynamic  loading  and 
structural  properties  of  such  faired  cables.  Additionally,  such  a study 
might  lead  to  an  optimized  design  for  future  faired  cable  systems  (optimized 
with  respect  to  load  carrying  characteristics  and  predictable  dynamics). 

In  this  report  we  will  document  the  development  of  a computer  program 
designed  to  study  precisely  those  hydrodynamic/structural  interactions  which 
contribute  to  the  "kiting"  of  faired  towlines.  Briefly,  the  present  effort 
has  been  directed  to  coupling  a collection  of  algorithms  which  predict  the 
hydrodynamic  loads  --  on  a fairing  cross-section  of  given  shape  --  to  a 
collection  of  algorithms  which  describe  the  fairing's  (structural)  response 
to  the  imposed  loading.  Naturally,  since  the  fairing  responds  to  loads  by 
deforming,  and  since  the  fairing's  shape  in  turn  determines  the  hydrodynamic 
loads,  an  equilibrium  state  is  achieved  only  after  several  "iterations" 
through  this  process.  Since  one  segment  of  the  program  computes  the 
HYDRQdynamic  loads  whose  effects  on  the  fairing  are  determined  by  the  non- 
linear Structural  Analysis  Program  segment,  we  have  (not  without  certain  mis- 
givings) called  the  full  program  HYDROSAP. 

The  report  is  organized  as  follows:  In  section  2 we  discuss  the  over- 
all organization  of  HYDROSAP.  Particular  attention  is  given  to  describing 
the  interface  and  auxiliary  software  necessary  to  effectively  control  all  the 


operations  of  the  program's  hydrodynamic  and  structural  response  modules. 

The  program's  input  requirements  are  explained  in  section  2.4.  In  section 
3 we  discuss,  in  detail,  the  hydrodynamic  computations  and  the  attendant 
(optional)  input  requirements.  A parallel  discussion  for  the  structural 
response  computations  appears  in  section  4.  Section  5 is  devoted  to  the 
description  and  discussion  of  a sample  HYDROSAP  analysis.  In  the  several 
appendices,  attached  hereto,  we  give  details  of  the  program's  structure  and 
describe  the  job  control  language  (JCL)  needed  to  operate  the  program  on  the 
IBM  360/91  system  at  the  Goddard  Space  Flight  Center. 


2.  DESCRIPTION  OF  THE  HYDROSAP  SYSTEM 


2.1  Overview 


The  HYDROSAP  program  has  been  built  to  study  the  nonlinear  interaction 
of  hydrodynamic  and  structural  effects  on  an  elastic  two-dimensional  fairing 
(cross-section)  exposed  to  a steady  field  of  fluid  flow.  In  order  to  achieve 
this  objective  elements  from  two  existing  computer  programs,  the  NASA  Multi  - 
component  Two-dimensional  Viscous  Airfoil  Program  (herein  dubbed  HYDRO)  and 
the  University  of  California  Nonlinear  Structural  Analysis  Program  (NONSAP), 
were  integrated  into  a single  software  system. 

The  integrated  system  determines  an  equilibrium  shape  for  a flexible 
fairing  exposed  to  the  hydrodynamic  loading  developed  by  moving  the  fairing 
through  a steady  state  stream.  The  equilibrium  shape  is  determined  by  pass- 
ing successively  through  the  HYDRO  and  NONSAP  modules  until  the  change  in 
fairing  geometry,  between  successive  passes,  is  smaller  than  some  pre-set 
error  tolerance.  This  iterative  computational  procedure  is  depicted 
pictorially  in  Figure  2-1. 

In  order  to  control  the  transfer  of  data  from  one  program  module  to  the 
other,  a collection  of  auxiliary  software  was  designed  and  coded.  These 
auxiliary  routines  perform  pre-  and  post-processing  operations  for  both 
modules.  They  also  provide  for  the  interfacing  between  modules.  In  addi- 
tion to  the  above,  the  HYDROSAP  package  includes  several  features  which  make 
the  program  quite  versatile,  but  yet  easily  accessible  to  (even)  the 
inexperienced  user: 

(1)  HYOROSAP  has  three  basic  modes  of  operation.  Besides  operating  in 
the  "full"  HYDROSAP  iteration  mode,  the  user  is  able  to  run  the 
HYDRO  or  NONSAP  modules  individually.  By  judiciously  overlaying 
(see  Appendix  A)  each  of  the  program  modules,  and  by  a careful 
design  of  the  auxiliary  interface  routines,  this  capability  has 
been  retained  with  little  increase  in  core  requirements  over  a 
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Figure  2-1.  Schematic  depicting  the  iterative  operations  of 

HYDROSAP  --  a hydrodynamic  loading,  structural  response 
software  system. 
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reduced  version  of  HYDROSAP  using  a "gutted"  form  of  the  HYDRO  and 
NONSAP  modules. 

(2)  The  various  HYDROSAP  options  and  features  are  activated  by  means 

of  a simple  NAMELIST  input.  This  form  of  input  frees  the  user  from 
the  error-prone  and  tiresome  preparation  of  formatted  input  decks. 
Moreover,  it  provides  a simple  and  effective  way  to  set  default 
values  for  control  variables. 

(3)  An  especially  attractive  HYDROSAP  feature  is  the  "short  form"  input 
pre-processor.  Both  the  HYDRO  program  and  the  NONSAP  program 
require  rather  large  and  complicated  formatted  input  decks.  If 
the  user  selects  the  built-in  short  form  input  pre-processor,  these 
input  decks  are  constructed  automatically,  based  on  the  values  of 
several  namelist  control  variables,  and  on  certain  assumptions 
about  the  class  of  symmetric  fairing  cross-sections  to  be  studied. 

Of  course,  the  user  may  elect  to  provide  complete  HYDRO  and  NONSAP 
decks  himself.  In  this  case,  the  fairing's  cross-section  (geo- 
metry) and  elastic  and  construction  details  can  be  completely 
arbitrary.  Whatever  the  origins  of  the  input  decks,  however,  the 
user  may  produce  an  "echo"  of  the  input  deck  card  images  on  the  out- 
put unit.  This  feature  facilitates  the  "debugging"  and  checking 

of  input. 

(4)  Also,  the  user  may  obtain  a summary  of  the  job  flow  at  his  remote 
terminal.  "Milestone"  messages  mark  the  beginning  and  ending  of 
each  stage  of  the  computation;  and,  convergence  messages  give  the 
"size"  of  the  displacement  error  (from  equilibrium)  relative  to  a 
user-assigned  convergence  tolerance.  Also,  any  program-generated 
error  messages  appear  at  the  terminal,  as  output,  in  addition  to 
appearing  on  the  main  printer  output. 

(5)  A convenient  checkpoint/restart  capability  is  also  available  to 
HYDROSAP.  Should  a run  terminate  abnormally  (e.o.,  the  computa- 
tion time  might  exceed  the  time  estimate  on  the  job  card)  the  user 
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may  restart  and  continue  the  computation  from  the  last  stored  check- 
point. This  procedure  represents  a savings,  in  time  and  costs, 
since  the  sequence  of  expensive  computations  completed  before  cut- 
off need  not  be  repeated. 

(6)  Finally,  data  necessary  to  produce  graphs  of  the  pressure  distribu- 
tion or  the  deformed  foil's  shape,  after  each  HYDROSAP  iteration, 
may  be  directed  to  disc  or  tape  storage  for  subsequent  retrieval 
by  the  BTS-developed  plotting  program  HSGRAF. 

2.2  HYDROSAP  Organization 


In  this  section,  we  describe  the  organization  of  the  HYDROSAP  system  and 
the  role  of  the  auxiliary  software  (Figure  2-2). 

The  HYDROSAP  computation  begins  (BEGIN0)  with  the  setting  of  default 
values  for  the  control  variables  in  NAMELIST  /DRV/.  INPUT0  reads  Tape  5 to 
ascertain  over-ride  control  variable  values  prescribed  for  the  run  by  the 
user.  A detailed  explanation  of  the  control  variables  is  given  in  section 
2.4.2  below. 

If  the  user  has  selected  short  form  input  pre-processing,  subroutines 
HYDIN  and  SAPIN  are  called  to  construct,  respectively,  the  HYDRO  (Tape  55) 
and  NONSAP  (Tape  66)  input  files.  The  fairing  cross-section  (profile)  is 
supplied  either  from  a collection  of  profile  library  subroutines  (FOIL01, 
FOIL02, • • -etc. ) by  DEFALT,  or  from  the  input  stream  through  NAMELIST 
/FAIRNG/.  Note  that  only  the  fairing's  upper  surface  profile  need  be  pro- 
vided, since  HYDIN  and  SAPIN  assume  a symmetric  fairing.  Control  variables 
from  NAMELIST  /DRV/  define  the  problem's  operating  conditions  (angle  of 
attack,  static  pressure,  and  freestream  velocity). 

Next,  SAPIN  generates  the  finite  element  mesh  for  NONSAP  based  on  the 
upper  surface  profile  geometry  and  on  values  of  the  control  variables  which 
define  the  mesh  "fineness",  and  the  materials  composition  for  the  fairing. 


START 


Figure  2-2.  HYDROSAP  ORGANIZATION  (and 

Computational  Path-Flow  Diagram) 


In  using  the  short  form  input,  the  fairing  is  assumed  to  consist  of  one  or 
more  of  the  following  parts:  a cable  core  (strength  member),  a fairing  tail 
(afterbody),  and/or  an  enclosing  skin.  Figure  2-3  shows  an  example  of  a 
finite  element  mesh  as  generated  by  SAPIN.  TAINT  is  a table  lookup  and 
interpolation  routine  used  here  by  SAPIN  to  convert  profile  surface  coordinates 
to  finite  element  surface  nodal  points. 

Of  course  the  user  may  decline  to  use  the  short  form  input  option,  in 
which  case  he  must  supply  the  complete  HYDRO  and  NONSAP  input  decks  on  Tape 
5.  (Refer  to  section  3.3  for  a detailed  description  of  the  HYDRO  input,  and 
to  section  4.3  for  a detailed  description  of  the  NONSAP  input.)  These  input 
decks  are  next  deconcatenated  and  stored  on  the  HYDRO  and  NONSAP  input  units, 

55  and  66,  respectively. 

With  pre-processing  completed,  HYDRO  is  called  to  compute  the  hydro- 
dynamic  loads  for  the  initial  fairing  shape.  In  addition  to  the  pressure 
coefficient  distribution,  HYDRO  provides  force  and  moment  coefficients  for 
the  profile.  Next,  CPTOP  converts  the  Cp  distribution  to  a pressure  load- 
ing distribution.  Subsequently,  TAINT  is  used  on  all  iterations  but  the 
first  to  update  the  fairing  profile  based  upon  the  deflections  computed  by 
NONSAP.  (The  HYDRO  module  is  discussed  at  length  in  section  3.) 

NONSAP  is  called,  in  turn,  to  compute  the  structural  deformation  induced 
by  the  pressure  loading.  Subroutines  LOADER,  NDLLDS,  and  TAINT  are  used  to 
convert  the  hydrodynamic  pressure  loading  to  an  (equivalent)  consistent 
finite  element  "nodal -load"  system.  As  explained  in  section  4,  NONSAP  itself 
employs  an  iterative  incremental  solution  procedure.  (We  refer  to  these 
NONSAP  iterations  as  "inner  loops".  The  iterations  between  HYDRO  and  NONSAP 
are  called  "outer  loops".) 

The  deformed  shape  is  next  passed  to  HYDRO,  which  computes  another  pres- 
sure distribution  based  upon  the  new  shape.  If  the  user  desires,  data  for 
graphics  post-processing  can  be  stored  on  Tape  88.  Another  program,  HSGRAF, 
is  available  to  produce  graphs  of  the  pressure  or  deformed  shape  profiles. 
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Cable  Fairing  Cross-Section  Finite  Element  Mesh,  as  Generated 
by  HYDROSAP  Short  Form  Input  Pre-processor 
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At  the  end  of  each  outer  loop,  CONVRG  is  called  to  determine  whether  an 
equilibrium  shape/pressure  distribution  has  been  obtained.  If  equilibrium 
is  not  achieved,  another  outer  loop  iteration  is  initiated.  Once  an  equilib- 
rium state  is  achieved,  or  a user-set  limit  on  the  number  of  outer  loop  itera- 
tions has  been  exceeded,  HYDROSAP  stops. 

Finally,  we  note  that  either  the  HYDRO  module  or  the  NONSAP  module  can 
be  completely  bypassed  should  the  user  wish  to  uncouple  the  hydrodynamic  and 
structural  computations. 

2 . 3 Interface  and  Auxiliary  Software 

In  this  section  we  will  give  short  descriptions  of  the  various  auxiliary 
and  interface  routines  which  are  used  in  the  HYDROSAP  system.  Complete  sub- 
routine descriptions  are  provided  in  Appendix  B. 

A subroutine  to  initialize  the  various  HYDROSAP  variables  and 
set  default  values  for  the  control  variables  in  NAMELIST  /DRV/. 

An  algorithm  to  test  the  size  of  the  current  relative  dis- 
placement increment  against  a user-assigned  tolerance.  This 
routine  is  called  after  the  NONSAP  displacement  computations 
are  completed,  within  each  HYDROSAP  outer  iteration. 

A subroutine  which  calls  the  appropriate  fairing  profile 
library  routine  during  short  form  input  pre-processing. 

The  main  HYDROSAP  executive  routine.  It  controls  execution 
of  the  HYDRO  and  NONSAP  program  modules  as  well  as  the  initial- 
ization procedure,  input  file  preparation,  convergence  testing, 
and  graphics  post-processing. 

This  is  a fairing  cross-section  upper  profile  library  routine 
numbered  "nm".  Up  to  99  such  library  routines  are  possible 
within  HYDROSAP.  These  routines  provide  default  surface  co- 
ordinates for  short  form  input.  Currently  nm  ranges  from  01 
to  05  as  follows: 
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FOILnm 


FOIL01 

FOIL02 

FOIL03 

FOIL04 

FOIL05 

HYDIN 


INPUT0 


OUTPST 


SAPIN 


TAINT 


LOADER 


NDLLDS 


CPTOP 


WRT14 


A typical  fairing  profile  measured  from  a sample  of  a DTNSRDC 
cable  section. 

Fairing  Model  I taken  from  DTNSRDC  Report  4610. 

An  NACA  0020  profile. 

An  NCS  0010  profile. 

A subroutine  reserved  for  future  expansion. 

This  routine  builds  a complete  HYDRO  input  file  based  upon 
values  of  the  control  variables  read  from  NAMELIST  /DRV/  dur- 
ing short  form  input  pre-processing. 

A routine  which  over-rides  the  default  values  of  the  HYDROSAP 
control  variables  with  values  read  from  NAMELIST  /DRV/. 

This  subroutine  writes  the  pressure  distribution  and  deformed 
shape  data  to  output  Tape  88. 

The  routine  which  builds  a complete  NONSAP  input  file  based 
upon  values  of  the  control  variables  from  NAMELIST  /DRV/  dur- 
ing short  form  input  pre-processing. 

This  is  a general  purpose  "table  lookup  and  interpolation" 
algorithm. 

An  algorithm  to  control  the  computation  of  consistent  nodal 
loads  which  are  to  be  applied  at  the  NONSAP  fairing  surface 
nodes,  based  on  the  computed  pressure  distribution  from  HYDRO. 

This  subroutine  determines  the  consistent  nodal  loads  using  the 
fairing's  surface  geometry  and  pressure  distribution. 

An  algorithm  which  converts  non-dimensional  surface  pressure 
coefficients  to  surface  pressures. 

This  subroutine  writes  the  HYDROSAP  checkpoint/restart  tape. 
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Tape  14.  WRT14  uses  TPSCAN  to  scan  the  four  NONSAP  restart 
tapes  to  determine  the  number  and  size  (in  bytes)  of  each 
tape's  records.  HYDRO  restart  data  is  concatenated  with  the 
four  NAPSAP  restart  tapes  to  form  each  HYDROSAP  checkpoint  on 
a single  unit.  Tape  14. 

TPSCAN  This  subroutine  scans  a file  written  with  unformatted  FORTRAN 

WRITE  statements  to  determine  the  length  (in  bytes)  of  each 
record  and  the  total  number  of  records. 

READ14  A routine  which  reads  and  unloads  the  checkpoint/restart  tape. 

Tape  14. 

2.4  HYDROSAP  Input 
2.4.1  Input  Stream  Format 

The  HYDROSAP  input  cards  are  identified  with  Tape  5 (usually  the  system 
input  stream).  The  input  stream  consists  of  six  distinct  card  groups,  not 
all  of  which  are  present  in  a single  run.  Only  the  first  two  card  groups 
must  be  present;  the  presence  or  absence  of  the  other  four  groups  is  deter- 
mined from  values  assigned  by  the  user  to  the  NONSAP  control  variables 
SHRTIN,  HYDSIM,  SAPS1M,  and  N.  All  of  the  control  variables  are  defined  in 
section  2.4.2,  below.  The  input  stream  card  groups  are  as  follows: 


CARD  GROUP  1 Title  Card 

Card  1 .1  Title  of  Run 

T 

Up  to  an  80  character  title  for  the  run. 
(This  card  must  appear.) 


CARO  GROUP  2 Control  Variables 

Card  2.1  $DRV 

12 

Header  card  for  NAMELIST  /DRV/;  signals  the  beginning  of 
the  control  input.  (This  card  must  appear.) 
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Card  2.2 


Card  2.n 
Card  2.n+l 


CARD  GROUP  3 


Card  3.1 


Card  3.2 


Card  3.m 


Card  5.m+1 


SHRTIN=T, • • • 

123 

Input  as  many  cards  as  required  to  override  default  control 
variable  values. 

SEND 

n 

End  card  for  NAMELIST  /DRV/;  signals  the  termination  of 
the  namelist  string.  (This  card  must  appear.) 

Cable  Fairing  Profile  Namelist  (optional) 

Card  Group  3 is  present  only  if  SHRTIN=T  and  N>0  in 
NAMELIST  /DRV/. 

SFAIRNG 

12 

Header  card  for  NAMELIST  /FAIRNG/. 

_ XX ( 1 )= • • • ,XX(N)= • • • , 

123 

ZZ ( 1 )= — ,ZZ(N)=- • • , 

123 

{(XX(I).ZZ(I)),  1=1, •••N)  are  the  profile  coordinates,  in 
inches,  for  the  upper  surface  of  the  symmetric  cable  fair- 
ing to  be  modeled.  N is  set  in  NAMELIST  /DRV/. 

SEND 

12 

End  card  for  NAMELIST  /FAIRNG/. 


CARD  GROUP  4 


Card  4.1 


HYDRO  Input  Deck  (optional) 

Card  Group  4 is  present  only  if  SHRTIN=F  and  SAPSIM=F  in 
NAMELIST  /DRV/. 

First  card  of  HYDRO  input  deck 
• } see  section  3.3 


Card  4.p  Last  card  of  HYDRO  input  deck. 
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CARD  GROUP  5 Spacer 

Card  5.1  99999999 

12345678 

A spacer  card  to  be  inserted  between  HYDRO  and  NONSAP  in- 
put decks.  (This  card  must  be  present  whenever  SHRTIN=F.) 


CARD  GROUP  6 


Card  6.1 

Card  6.q 


NONSAP  Input  Deck  (optional) 

Card  Group  6 is  present  only  if  SHRTIN=F  and  HYDSIM=F  in 
NAMELIST  /DRV/. 

First  card  of  the  NONSAP  input  deck. 

• } see  section  4.3 
Last  card  of  the  NONSAP  input  deck. 


2.4.2  Control  Variable  Definitions 


The  HYDROSAP  Control  Variables  (Card  Group  2 of  the  HYDROSAP  Input)  fall 
into  two  groups:  logical  variables  and  arithmetic  variables.  The  logical 
variables  are  true/false  switches  used  to  activate  the  various  program  options. 
The  arithmetic  variables  are  used  to  set  values  for  the  integer  and  real 
variables  which  control  the  flow  of  computations  (examples  are:  iteration 
counters,  output  units,  mesh  fineness,  and  convergence  tolerance),  or  define 
the  physical  parameters  of  the  problem  under  study  (such  as  flow  conditions 
and  material  properties).  Detail  definitions  of  the  control  variables  are 
as  follows: 


VARIABLE 

CHKPNT 

(CHeckPoiNT) 

CORE 


TYPE  DEFAULT  INTERPRETATION 


LOGICAL  .FALSE.  =.TRUE.  Produce  checkpoint 

tape  (FT14F001) 

=. FALSE.  Checkpointing  dis- 
abled. 

LOGICAL  .FALSE.  =.TRUE.  Fairing  has  a cable 

(core)  material 
=. FALSE.  No  distinct  cable 
(core)  present 
(SHRTIN= .TRUE . only) 
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VARIABLE 

TYPE 

DEFAULT 

INTERPRETATION 

DEBUG 

LOGICAL 

.FALSE. 

=.TRUE.  Special  print  in 

LOADER:  HYDRO  and 
NONSAP  output  pro- 
vided at  each  HS 
iteration  (itera- 
tion = outer  loop) 

=. FALSE.  HS  print  produced  at 
start  of  run  only. 

No  debug  output. 

ECHO 

LOGICAL 

.TRUE. 

= .TRUE . Echo  of  HYDRO  and 
NONSAP  input  decks 
produced. 

= . FALSE.  Suppress  ''echo"  of 
files . 

HYDSIM 

(HYDro  SIMulation) 

LOGICAL 

.FALSE. 

= .TRUE . Program  operates  in 
the  HYDRO  simulation 
mode  (SAPSIM  must  be 
.FALSE.) 

=. FALSE.  Signifies  a request 

for  the  full  HYDROSAP 
run;  or  NONSAP 
simulation,  depend- 
ing on  the  value 
assigned  to  SAPSIM. 

POST 

LOGICAL 

.FALSE. 

=.TRUE.  The  Graphics  Post- 
processor output  is 
directed  to  Unit  88. 

=. FALSE.  No  operation  requested 

RESTRT 

LOGICAL 

.FALSE. 

= .TRUE . Read  checkpoint  tape 
and  restart  at  time 
TSTART 

=. FALSE.  Restart  disabled. 

SAPS  I M 

LOGICAL 

.FALSE. 

=.TRUE.  Program  is  to  operate 
in  NONSAP  simulation 
mode  (HYDSIM  must  be 

set  to  .FALSE7T" 

=. FALSE.  Full  HYDROSAP  run; 

or  HYDRO  simulation, 
depending  upon  value 
assigned  to  HYDSIM. 


VARIABLE 

TYPE 

DEFAULT 

INTERPRETATION 

SHRTIN 

(SHoRT  form  INput) 

LOGICAL 

.TRUE. 

=.TRUE.  Signifies  the  short 
form  input.  An  auto- 
matic input  file  gen- 
eration occurs  (for 
foil  geometry  and 
state) . 

=. FALSE.  HYDRO  and  NONSAP  in- 
put decks,  separated 
bv  999999999,  must 
folT5w'$BR^-$END  on 

SYSIN  file. 

SKIN 

tt 

LOGICAL 

.FALSE. 

=.TRUE.  Denotes  the  fairing 
cross-section  has  a 
skin  attached. 

=. FALSE.  No  skin  member  is 
presumed . 

TAIL 

( 

1 

LOGICAL 

.TRUE. 

= . TRUE . The  fairing  has  a 

fairing  "tail"  mem- 
ber. 

=. FALSE.  No  fairing  tail  is 

assumed  (CORE  must  be 
set  to  .TRUE,  in  this 
case) . 

TPRINT 

(Terminal  PRINT) 

LOGICAL 

.TRUE. 

=.TRUC.  Terminal  summary  out- 
put (error  msgs.  and 
job  flow) 

=. FALSE.  No  terminal  print. 

INTORD 

(INTerpolation  ORDer) 

INTEGER 

3 

Denotes  the  interpolation 
order  for  the  pressure  dis- 
tribution conversion  from  the 

HYDRO  surface  to  the  NONSAP 
surface. 

INTORX 

(INTerpolation  ORder  X) 

INTEGER 

1 

Interpolation  order  to  be 
assigned  for  the  displacement 
conversion  from  the  NONSAP 
surface  to  the  HYDRO  surface. 

I PANEL 

f | 

INTEGER 

2 

IPANEL  is  equivalent  to  HYDRO 
control  variable  IOP. 

IPAPER 

INTEGER 

6 

Unit  number  assigned  to  the 

HYDROSAP  DRIVER  module  output. 
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VARIABLE 

TYPE 

DEFAULT 

INTERPRETATION 

IPRINT 

INTEGER 

=NSTEPS 

Denotes  the  NONSAP  print  fre- 
quency within  HYDROSAP 
ITERATION  #1. 

ISMTH 

INTEGER 

0 

ISMTH  is  equivalent  to  HYDRO 
control  variable  ISMO. 

ITRSWT 

(ITeRation  SWiTch) 

INTEGER 

8 

The  number  of  HYDROSAP  itera- 
tions using  HYDRO  potential 
flow  solution  mode  only. 

IWRT0 

INTEGER 

6 

Denotes  the  NONSAP  output 

UNIT. 

IWRT1 

INTEGER 

6 

Denotes  the  HYDRO  output  unit. 

LOOPMX 

(LOOPMaXimum) 

INTEGER 

1 

Signifies  the  maximum  number 
of  HYDROSAP  iterations  (outer 
loops)  allowed. 

MESHC 

(MESHChordwise) 

INTEGER 

1 

Denotes  the  number  of  chord- 
wise  interior  stations  used 
to  generate  the  finite  ele- 
ment mesh,  when  SHRTIN  is 
set  to  .TRUE. 

MESHT 

(MESHThickness-wise) 

INTEGER 

1 

Denotes  the  number  of  lateral 
interior  stations  used  to  qen- 
erate  the  finite  element  mesh 
when  SHRTIN  is  set  to  .TRUE. 

N 

INTEGER 

0 

N>0  Signifies  the  number  of 
upper  surface  coordinate 

pairs  to  be  input  through 
NAMELIST  $FAIRNG-$END,  fol- 
lowing NAMELIST  SDRV-SEND, 
on  SYSIN. 

N=0  Means  no  fairing  profile 
is  input. 

N<0  A pointer  for  profile 
1 ibrary  routine  |N| . 

Thus,  N=  -2  reads  profile 
points  from  subroutine 
FOIL02. 

(Applicable  when  SHRTIN  set  to 

.TRUE.,  only.) 


VARIABLE 

TYPE 

DEFAULT 

INTERPRETATION 

NSTEPS 

INTEGER 

20 

Denotes  the  number  of  (quasi- 
static) load  steps  to  be  used 
in  applying  the  field  pressure 
load  to  a fairing,  in  NONSAP, 
during  HYDROSAP  iteration  # 1. 

Note: 

NSTEPS 

LOADS 

should  be  increased  when  the  message  "OUT-OF-BALANCE 
EXCEED  INCREMENTAL  LOADS"  is  obtained. 

ALPHA 

REAL*8 

0.D0 

Angle  of  attack  (deg. ) 

ECORE 

REAL *8 

1.D06 

Young's  modulus  for  the  core 
material  (psi). 

ESKIN 

REAL *8 

1 .D04 

Young's  modulus  for  the  skin 
material  (psi). 

ETAIL 

REAL*8 

1 . D04 

Young's  modulus  for  the  fair- 
ing tail  material  (psi). 

PCORE 

PSKIN 

P TAIL 

REAL*8 

.3 

0 

.3 

Poisson  ratio  for  the  core, 
skin,  and  tail, 
respectively. 

PZERO 

REAL*8 

14.7D0 

Static  pressure  loading  (psi). 

RHO 

REAL*8 

2.D0 

H^O  density  (slugs/ft3). 

TOL 

REAL*8 

1*10'2 

Denotes  the  HYDROSAP  conver- 
gence tolerance  level.  Con- 
vergence is  achieved  when  rms 
displacement  increment 
<(rms  displacement  ) * TOL. 

TSKIN 

REAL*8 

.005 

Skin  thickness  (inches) 

(applies  only  if  SKIN  is  set 
to  .TRUE.). 

TSTART 

(Time  of  reSTART) 

LOGICAL 

0.D0 

Applies  only  when  RESTRT= 

.TRUE.  HYDROSAP  will  restart 
at  checkpoint  for  time  TSTART. 

VFRSTR 

REAL*8 

22.  D0 

Free  stream  velocity  (fps). 

Y1CORE 

REAL *8 

0.D0 

Chordwise  location  of  the  core 
starting  point. 

Y2CORE 

REAL*8 

0.D0 

Chordwise  location  of  the  core 
ending  point. 
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2.5  HYDROSAP  Output 

The  output  for  the  entire  HYDROSAP  program  can  be  divided  into  four 
groups:  HYDROSAP  printer  output,  HYDROSAP  terminal  summary  output,  HYDRO 
module  output,  and  NONSAP  module  output.  The  user  may  control  the  presence 
of  these  output  forms  through  the  HYDROSAP  control  variables  IPAPER,  IWRT0, 
IWRT1 , and  TPR1NT,  as  discussed  in  section  3.4  and  section  4.4  below.  In 
this  section  we  will  briefly  discuss  the  HYDROSAP  printer  and  terminal  out- 
put. 

2.5.1  HYDROSAP  Printer  Output 

HYDROSAP  printer  output  appears  on  unit  IPAPER  (usually  unit  6). 

Basically  it  contains  a listing  of  program  control  information,  convergence 
data  regarding  the  HYDRO-NONSAP  iterations,  for  "full"  HYDROSAP  runs,  and 
any  error  messages  which  may  occur. 

Immediately  following  the  banner  title  is  an  echo  of  the  values  of  the 
entire  control  variable  namelist,  /DRV/.  If  the  user  has  set  ECH0=T , the 
next  data  appearing  on  the  printer  listing  is  an  echo  of  the  HYDRO  and/or 
NONSAP  input  deck  card  images.  The  echo  is  preceded  by  the  message 

ECHO  OF  INPUT  DECK  FOR  HYDRO  AND  NONSAP  SUBPROGRAMS 
and  followed  by  the  message 

END  ECHO  . 

The  HYDRO  and  NONSAP  decks  are  separated  by  the  notation 
--END  HYDRO  INPUT;  BEGIN  NONSAP  INPUT— 

Following  the  input  stream  echo  on  unit  IPAPER  are  messages  marking  the 
beginning  of  each  outer  iteration.  The  usual  practice  is  to  direct  HYDRO  and 
NONSAP  output  to  the  printer  as  well  ( I WRT1  = I WRT 0=IPAPER=6) , thus  the  results 
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of  the  hydrodynamic  and  structural  deformation  computations  appear  in  sequence 
for  each  HYDROSAP  outer  iteration. 

Following  the  HYDRO  and  NONSAP  output  for  each  iteration  is  a message 
from  subroutine  CONVRG.  The  current  relative  displacemei.t  error,  RATIO,  and 
the  user's  input  convergence  tolerance,  TOL,  are  listed.  ISTOP  equals 
zero  until  RATIO^TOL,  in  which  case  ISTOP  is  set  to  1 and  the  computations 
are  terminated.  The  computations  also  terminate  if  the  iteration  counter, 
ITER,  exceeds  the  user-set  maximum,  LOOPMX.  The  displacement  increment  and 
total  displacement  for  each  degree  of  freedom  are  listed  in  compacted  form 
in  order  of  equation  number.  Finally,  the  output  ends  with  an  informative 
message  explaining  how  the  computation  terminated. 

2.5.2  HYDRQSAP  Terminal  Summary  Output 

Terminal  summary  output  appears  on  unit  12  if  TPRINT=T.  A sample  of 
terminal  summary  output  is  given  in  Figure  2-4.  The  user  is  informed  as  the 
HYDRO  and  NONSAP  modules  are  accessed  during  each  HYDROSAP  outer  iteration. 

The  results  of  the  convergence  check,  at  each  outer  iteration,  are  also  given; 
and,  a final  message  informs  the  user  of  the  run's  method  of  termination. 

Under  abnormal  circumstances  any  error  messages  generated  by  the  program  are 
reproduced  here  also. 


i 
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,4 

X 


!' 

- 

■COMMENCE  HYDPOSRP  EXECUTION 
♦♦♦  ITERATION  1 ♦♦♦ 

BEGIN  HYDPD  PHASE 
END  HYDPO  PHASE 
BEGIN  NDNSRP  PHRSE 
END  NONSRP  PHRSE 

ERROR  = I.OOOOD  00.  IJ.P.T.  UNITY  = 1 . 0000D  03*  NOT  CONVERGED 

♦♦♦  ITEPRTION  £ 

BEblN  HYDRO  PHRSE 
END  HYDPO  PHRSE 
BEGIN  NONSRP  PHRSE 
END  NONSRP  PHRSE 

EPPOP  = 8.5004D-05*  W.R.T.  UNITY  = 8.5004D-08.  CONVERGED 

PUN  TERM1NRTED  BECRUSE  CONVERGENCE  ACHIEVED 


Figure  2-4.  Sample  of  Terminal  Print  Summary 
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3.  DESCRIPTION  OF  HYDRO  - THE  VISCOUS/POTENTIAL  FLOW  SUBPROGRAM 

This  component  of  HYDROSAP  is  the  program  segment  used  to  describe 
the  hydrodynamic  loads  which  are  applied  to  a streamlined  fairing's  cross- 
section.  The  particular  collection  of  algorithms  and  subroutines,  found 
here,  have  evolved  from  the  systems  described  in  [1,2],*  and  provide  a basis 
for  the  system  described  in  [2].  The  utility  of  this  program  segment  has 
been  alluded  to  earlier;  therefore,  the  present  section  of  the  documentation 
has  been  prepared  to  provide  the  reader  with  a general  understanding  of  the 
program's  operation,  its  attributes  and  its  general  limitations. 

Included  in  the  following  paragraphs  is  an  overview  and  description  of 
HYDRO;  a short  description  of  the  subroutines,  etc.  comprising  the  system; 
and  a brief  discussion  on  the  input  format  required  to  run  the  program. 

Should  the  reader  desire  a more  detailed  discussion  and  description  of  the 
theory,  etc.,  used  in  building  the  system,  the  references  (included  herein) 
should  be  consulted  (especially  [1],  [2],  [3]  and  [4]). 

3.1  Program  Description 

The  "hydrodynamic  loading"  program,  which  has  been  given  the  name  HYDRO, 
here,  is  constructed  in  four  segments.  Each  of  these  subsystems  plays  a dis- 
tinct computational  role  in  developing  the  foil's  pressure  distribution, 
and  in  the  generation  of  force  and  moment  coefficients  which  characterize 
the  fairing  section.  A descriptive  flow  chart  for  this  program  is  given 
on  Figure  3.1 . 

The  principal  program  segments  here  are  (each)  controlled  by  one  of  the 
following  subroutines:  AIRGM,  POTFL,  BDLYR  and  LOADS0.  In  the  next  few 
paragraphs  a word  description  of  their  operation  and  function  will  be  given. 
Needless  to  say,  there  is  one  overriding  (executive  ) routine,  in  the  pro- 
gram's hierarchy,  which  has  the  function  of  initiating,  maintaining  and 
executing  control  over  the  entire  subsystem.  This  routine  has  been  (appro- 
priately) dubbed,  HYDRO. 


START' 


INPUT 

READ:  control  information  and 
input  data 

PHASE 

Compute  geometry  for  panel  distribution. 

AND 

loft  foil  and  determine  longest  chord 

CONTROL 

OPERATION 

Establish  e-distribution  for  foil 

geometry 

POTENTIAL 

FLOW 

(FIELD) 

COMPUTATIONS 


Compute  velocity  distribution  due 
to  foil  thickness 

Compute  velocity  distribution  due 
to  camber  (only) 


VISCOUS 

FLOW 

(FIELD) 

COMPUTATIONS 


♦ 

i 

“ “1 

IF  IPOT  = (0,1)  viscous  boundary  layer  is 

simulated  by  (geometry  change,  distributed 

| 

source  singularities) 

1 

Combine  velocities,  (and  for  subsequent 

| 

computational  manipulations)  describe  viscous 

velocity  distribution  by  iterative  operations 

| 

(based  on  boundary  layer  influences) 

1 

1 

PRESSURE 

DISTRIBUTION 

DETERMINATION 


FORCE  AND 

MOMENT  COEFFICIENTS, 
AND  PRESSURE 
LOADING 
COMPUTED 


Compute  Pressure  Distribution  for 
foil  geometry  (and  conditions) 
assumed 


Compute  hydrodynamic  Force  and  Moment 
coefficients  (from  known  pressure  and 
shear  (force)  distributions) 

Calculate  pressure  loading,  for  the 
foil,  based  on  known  pressure 
distribution  and  flow  field  state 


Figure  3.1.  HYDRO  Flowchart.  Program  Operation  for 
Single  Element  Streamlined  Foil 
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3.1.1  The  Foil  Geometry  Section 


The  AIRGM  segment  is  responsible  for  reading  in  the  input  data,  perform- 
ing the  various  geometry  chores  needed  in  setting  up  a problem,  and  establish- 
ing the  various  geometric  relationships  required  in  the  subsequent  analysis. 
One  of  the  major  chores  performed  by  this  routine  is  the  distribution  of 
(geometry)  "panels"  needed  for  the  pressure  distribution  calculations.  In 
this  regard  the  system  has  three  options  available  internally;  the  choice 
of  option  is  user  controlled  through  the  input  stream.  The  default  option 
(presently  available)  is  a distribution  of  "panels"  according  to  the  degree 
of  local  surface  curvature. 

The  various  subroutines  making  up  the  AIRGM  segment  of  this  program  can 
be  seen  on  the  HYDRO  structure  map  in  Appendix  A. 

3.1.2  The  Potential  Flow  Field  Section 


The  next  segment  called  is  that  which  computes  the  potential  flow  solu- 
tions for  the  foil  geometries  prescribed.  The  primary  subroutine  here  is 
POTFL;  the  other  subroutines  comprising  this  segment  are  readily  identified 
on  the  structure  map  (Appendix  A). 

In  the  calculations  performed  here  the  velocity  distribution,  for  a foil, 
is  developed  by  a superposition  of  velocity  fields.  In  this  regard  the 
velocity  contribution  due  to  foil  thickness  (alone)  is  added  to  that  due  to 
camber,  and  this  combination  is  subsequently  modified  to  show  the  influence 
of  boundary  layer  (as  it  develops  over  the  surfaces).  In  pursuit  of  this 
latter  effect, the  program,  controlled  by  input,  has  the  option  of  simulating 
the  boundary  layer's  effects  by:  (1)  modifying  the  foil's  geometry;  or, 

(2)  introducing  a distribution  of  source  singularities  to  account  for  the 
layer's  buildup. 

During  the  first  two  computational  iterations  (in  this  program  segment) 
the  velocity  distributions  due  to  thickness  and  camber  --  for  the  basic  foil 
geometry  --  are  computed.  For  all  subsequent  iterates  these  velocities  are 
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recalculated,  accounting  for  the  influence  of  boundary  layer  displacement 
thickness  over  the  basic  foil's  geometry.  Of  course,  when  the  viscous 
boundary  layer  is  neglected  the  computational  procedure  is  considerably 
shortened  and  the  need  for  "iterations"  is  removed. 

Calculations  for  the  velocity  fields  (due  to  foil  thickness  and  camber) 
make  use  of  Oellers'  Vortex  Distribution  Method  [5],  using  the  scheme  which 
is  outlined  in  [4].  Incorporated  in  with  these  computations  is  the  procedure 
used  (here)  to  satisfy  the  Kutta  condition  for  the  foils.  The  procedure 
used  represents  a modified  Kutta  condition  wherein  the  upper  and  lower  sur- 
face vortex  strengths  are  matched  in  magnitude,  but  assigned  opposing 
signs. 

If  the  viscous  boundary  layer  is  simulated  by  a distribution  of  source 
singularities  then  the  subroutine  (SOURCE)  is  called.  Otherwise  the  bound- 
ary layer's  influence  is  provided  for  by  modifying  the  velocity  field  for 
foil  camber.  For  those  cases  where  trailing  edge  thickness  is  not  zero,  the 
velocity  field  due  to  thickness  is  altered  through  a closure  of  the  geometry 
one  chord  length  downstream  of  the  true  trailing  edge. 

One  difficulty  with  the  present  program  is  that  it  does  not  have  the 
ability  to  handle  separated  flow  (until  recently*,  this  was  typical  of  the 
various  analytical  simulations).  Another  discrepancy,  which  may  be  remedied 
soon,  is  the  manner  in  which  the  wake  (behind  the  foil)  is  modeled;  and, 
the  attendant  drag  increment  attributed  to  such  a wake.  The  present  system 
accounts  for  the  increase  in  drag  solely  through  an  empirical  procedure  based 
on  tests  with  truncated  airfoils. 

Since  the  velocity  field,  about  the  foil,  is  determined  here  then  the 
pressure  distribution  is  computed  here,  also.  In  addition,  should  there  be 
any  reason  to  make  corrections  for  compressibility,  then  this  would  be  the 
place  to  do  so.  (For  those  readers  desiring  a concise  but  illuminating  dis- 
cussion on  various  aspects  of  the  general  theories  alluded  to,  here,  ref- 
erence [6]  is  recommended.) 

*See  Reference  [15] . 
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3.1.3  The  Boundary  Layer  Computations  Section 


The  boundary  layer  segment  of  the  HYDRO  program  is  designed  to  provide 
needed  calculations  for  the  boundary  layer's  influence.  The  segment  is 
controlled  by  the  routine  BDLYR,  with  calls  being  made  (primarily)  to  LAMBl 
and  TURBL.  There  the  laminar  and  turbulent  boundary  layer  algorithms 
reside;  and  there  the  needed  caluclations  are  made.  Since  HYDROSAP  is 
designed  to  analyze  a single  element  foil  (only)  the  confluent  boundary 
layer  portions  of  this  segment  are  not  called  or  exercised. 

The  subroutine  LAMBL  is  based  on  the  theories  of  Cohen  and  Reshotko 
[7],  as  modified  by  Goradia,  and  work  published  in  Schlichting  [8].  In 
addition  to  computing  the  boundary  layer  thickness,  the  program  determines 
the  (laminar)  friction  coefficients  and  establishes  laminar  separation  loci. 

The  subroutine  BLTRAN  is  called  to  determine  the  boundary  layer  trans- 
ition points  (if  the  free  transition  option  is  prescribed). 

The  turbulent  boundary  layer  calculations  are  provided  for  in  the  sub- 
routine TURBL.  There,  the  algorithm  (used)  is  based  on  Truckenbrodt ' s 
Integral  Method  [9],  as  modified  by  Goradia.  In  addition  to  computing  the 
usual  viscous  flow  parameters,  this  subroutine  also  calculates  a value  for 
the  drag  coefficient  based  on  the  approximate  theory  due  to  Squire  and 
Young  [10].  (See  the  structure  map.  Appendix  A,  for  details  regarding  the 
other  subroutines  contributing  to  the  boundary  layer  segment  of  the  pro- 
gram. ) 

3.1.4  The  Hydrodynamic  Loading  Parameters  Section 

The  last  segment  of  this  subprogram  is  controlled  by  the  subroutine 

LOADS0.  Here  the  (usual)  force  and  moment  coefficients  are  computed,  using 

the  known  pressure  and  shear  force  distributions  existing  on  the  foil.  In 

this  algorithm  the  incremental  coefficients  are  evaluated  for  contributions 

due  to  C and  due  to  C,  ; these  are  then  converted  to  C,  and  Cn 
p f L D 

values;  the  moment  coefficients  at  (near)  the  aerodynamic  center  (C/4)  are 
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determined;  and,  the  lift-to-drag  ratio  is  established. 

With  these  data  in  hand  the  engineer  is  made  aware  of  the  foil's 
sectional  hydrodynamic  characteristics  --  as  these  are  computed  for  the 
deflected  shape  (geometry)  --  and,  the  last  operational  step  to  be  performed 
is  the  conversion  of  these  loading  parameters  to  pressure  loads.  The  sub- 
routine CPTOP  is  called  for  this  final  conversion  operation.  After  this 
last  task  is  completed,  the  loads  can  be  passed  to  the  NONSAP  subprogram 
where  the  foil's  elastic  response  (deflections)  will  be  ascertained. 

3.2  Subroutine  Descriptions 

In  this  section  of  this  document  a brief  word  description  of  all  sub- 
routines comprising  the  HYDRO  Subprogram  is  presented.  This  information, 
coupled  with  the  chart  shown  in  Appendix  A,  will  add  to  the  statements 
above,  acting  to  provide  a general  understanding  of  how  the  subprogram  func- 
tions. 

AIRGK  Reads  input  data  for  a (multi-component)  airfoil,  and  does  the 
computations  needed  to  describe  the  geometry  of  the  Panel  Model 
used  in  the  Potential  Flow  and  Boundary  Layer  Programs. 

BDLYR  This  subprogram  serves  as  a control  routine  performing  the  manage- 
ment chores  for  all  Boundary  Layer  computations. 

BLTRAN  The  Boundary  Layer  Transition  routine;  performs  those  computations 
needed  to  determine  and  describe  boundary  layer  transition 
phenonema . 


CAMBER 


An  algorithm  which  computes  the  coordinates  for  the  "equivalent" 
airfoil;  an  airfoil  modified  to  account  for  boundary  layer  dis- 
placement thickness.  The  new  (equivalent)  airfoil  shape  is 
obtained  by  adding  a camber  line  (addition)  computed  by  averaging 
the  upper  and  lower  surface  displacement  thicknesses  to  the  (base 
camber  line  for  the  basic  airfoil  geometry. 


) 
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COMPR 


This  subprogram  computes  the  local  pressure  coefficient  (distribu- 
tion) using  the  Karman-Tsien  correction  (law);  also,  it  estab- 
lished the  stagnation  point  locus  (on  the  airfoil). 


CONFBL*  A subprogram  which  provides  the  computational  needs  in  calculat- 
ing the  confluent  boundary  layer. 

C0NF5*  Goradia's  computational  method  for  a confluent  boundary  layer, 
in  the  "core  region". 

C0NF7*  Goradia's  computational  method  for  the  confluent  boundary  layer, 
in  the  "main  region"  of  the  airfoil  system  (main  region  1). 

C0NF8*  Goradia's  computational  procedure  for  the  confluent  boundary 
layer  (in  main  region  2). 

CUBSPL*  A computational  algorithm  used  to  fit  a cubic  spline  through  a 
curve  form,  y(x)  . 

CURV  (A  Function  routine.)  Employed  to  compute  surface  curvature,  at 

a given  airfoil  surface  location  (at  a specific  Theta  locus). 

DERV  A subroutine  which  computes  derivaties  of  a tabulated  array, 

Y - F ( X ) . The  method  utilized  a weighted  average  of  derivatives 
obained  by  a quadratic  fit  to  points  (J-2,  J-l,  J),  (J-l,  J,  J+l ) 
and  (J,  J+l,  J+2). 

DLIM  (A  specialize  function  (routine)  called  by  CONFBL.) 

DSTEQ  A computational  procedure  used  to  "smooth"  the  DSTAR  distribution 
output  from  the  boundary  layer  routines. 


*These  subroutines  are  not  required  for  the  single  element  foils  dealt  with 
in  the  present  study;  however,  they  are  retained  (and  discussed)  since  they 
constitute  a part  of  the  total  computational  program  and  may  be  of  value 
to  future  users. 


4 
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PROOT  A subprogram  which  finds  the  roots  of  a cubic  polynomial.  In 
this  routine:  A = a coefficient  vector,  and  X = a solution 
vector,  for  A(1)+A(2)*X+A(3)*X**2+A(4)*X**3=0;  X is  complex 
and  X = X ( 1 )+lX(2) . 

PRTCRD  A "print"  routine  --  "prints"  the  coordinates  of  a lofted  multi- 
component  airfoil  system. 

PRTGM  A "print"  routine  --  "prints"  the  input  airfoil  geometry,  prior 
to  lofting. 

READIT  A "read"  routine  --  "reads"  the  input  data  for  the  multi-component, 
two-dimensional  viscous  airfoil  programs. 

S1MS0L  A "solution"  algorithm  --  solves  the  linear  simultaneous  set  of 
equations,  AX  = B,  in  N-unknowns,  and  returns  the  solution  in  B. 

SMOCR  A "smoothing"  routine  --  smooths  the  input  airfoil  coordinates  by 
consecutively  fitting  a least  squares  polynomial  (of  degree  four) 
through  seven  points  at  a time. 

SMOOTH  A subprogram  designed  to  smooth  an  array,  Y(X),  wherein  Y(X)  is 
replaced  by  a weighted  average  of  Y(X)  and  the  average  of  two  (2) 
predictions  of  the  Y(X)  values  obtained  from  a least  squares  fit 
to  four  (4)  values  of  Y(X)  beginning  at  the  jth  and  - 1 loci, 
respectively. 

SOURCE  A subroutine  developed  to  compute  the  source  singularities 

(strengths)  which  will  produce  the  equivalence  to  a boundary  layer 
displacement  thickness. 

THKVL  An  algorithm  designed  to  compute  a velocity  distribution  based  on 
airfoil  thickness  for  a viscous  flow  airfoil  shape  oriented  at 
zero  degress  angle  of  attack. 
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GAPA*  This  routine  computes  the  "gap"  distance  between  forward  and  after 
components  of  a mul ti -component  airfoil  system.  Gap  is  defined 
as  the  distance  from  a lower  surface  trailing  edge  to  the  upper 
surface  of  a following  (after  airfoil)  element. 

GEOM  This  routine  lofts  and  "panels"  each  component  of  a multicomponent 
airfoil  system. 

HYDRO  The  main  control  (management)  routine  for  the  entire  computational 
procedure.  This  subprogram  makes  calls  to:  AIRGM  (for  input 
data  and  geometry  calculations  pertaining  to  the  aerodynamics 
model);  POTFL  (where  the  potential  flow  calculations  are  executed, 
about  the  equivalent  airfoil  shape(s));  BDLYR  (where  boundary 
layer  properties  are  established);  LOADS0  (where  the  aerodynamic 
loads  are  computed);  and  to  the  various  graphics  routines  (where 
graphics  output  data  are  stored,  and  readied  for  presentation 
via  a selected  mode  of  presentation). 

INPUT**  A subprogram  which  lists  input  requirements  for  the  multi-component, 
two-dimensional  viscous  airfoil  program. 

INTER  An  interpolat ion  routine. 

KUTTA  An  algorithm  which  computes  the  aerodynamic  influence  coefficients 
for  a modified  Kutta  condition.  (The  condition,  as  stated, 
requires  that  the  unit  vortex  strength(s)  vary  quadratically  over 
the  last  four  (4)  panel  (end)  points  near  to  the  upper  and  lower 
surface  trailing  edge  loci;  and,  that  at  the  trailing  edge  the 
upper  and  lower  surface  unit  vortex  strengths  be  equal  in  magni- 
tude but  opposite  in  sign.) 


*These  subroutines  are  not  required  for  the  single  element  foils  dealt  with 
in  the  present  study;  however,  they  are  retained  (and  discussed)  since 
they  constitute  a part  of  the  total  computational  program  and  may  be  of 
value  to  future  users. 

**Not  present  in  the  current  version  of  HYDROSAP. 
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LAMBL 


The  laminar  boundary  layer  algorithm. 

LOADS0  This  subprogram  computes  the  force  and  moment  coefficients,  for 
the  mul ti -component  airfoils,  using  the  known  pressure  and  shear 
(force)  distributions. 

OELLER  Computes  the  aerodynamic  influence  coefficients,  for  the  airfoils' 
potential  flow  field,  using  Oeller's  singularity  distribution 
method  (see  NASA  CR-2523  (1975)  for  a detailed  description  of  the 
method). 

OVRLAP*  A computational  algorithm  which  calculates  the  overlap  distance 
between  "forward"  and  "aft"  components  of  a mul ti -component  air- 
foil system. 

PANEL  A subprogram  designed  to  distribute  a desired  number  of  (computa- 
tional) panels  over  the  (profile)  surface  of  an  airfoil  shape. 

POTFL  An  algorithm  for  computing  the  viscous/potential  flow  about  an 
equivalent  mul ti -component  airfoil  shape  --  one  composed  of  the 
basic  airfoil  (geometry)  with  the  boundary  layer  displacement 
thickness  added. 

Note:  If  IPOT  = 0,  the  viscous  boundary  layer  is  simulated  as 
change  in  geometry,  IPOT  = 1,  the  boundary  layer  effect 
is  simulated  by  distributed  source  singularities. 

POTVL  This  subroutine  calculates  a potential  flow  velocity  distribution, 
over  a multi -component  airfoil  system,  using  Oeller's  Vortex 
Distribution  Method,  but  with  boundary  conditions  modified  to 
include  a distributed  source  field. 


♦These  subroutines  are  not  required  for  the  single  element  foils  dealt  with 
in  the  present  study;  however,  they  are  retained  (and  discussed)  since 
they  constitute  a part  of  the  total  computational  program  and  may  be  of 
value  to  future  users. 
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TRANS*  A subroutine  which  is  employed  to  translate  and  rotate  the  co- 
ordinates of  the  several  airfoil  components,  of  a multi -component 
system,  relative  to  the  basic  frame  of  reference. 

TURBL  An  algorithm  developed  to  compute  the  turbulent  boundary  layer 
properties  via  the  Truckenbrodt  Integral  Method  (as  modified  by 
Goradia) . 

3.3  HYDRO  Input 

This  section  gives  details  on  the  required  input  file  content  to  the 
HYDRO  subprogram  module.  When  the  HYDROSAP  program  is  operating  with 
SHRTIN  = .TRUE.,  the  short  form  input  mode,  this  file  is  automatical ly 
developed  internally.  On  the  other  hand,  when  a specialized  user  applica- 
tion program  is  to  be  employed,  the  user  will  necessarily  prepare  (much  of) 
this  file  and  read  it  in  as  program  input. 

No  discussion,  or  description,  of  those  control  cards  which  are  not 
applicable  to  the  present  use  of  this  program  (e.g.,  the  confluent  boundary 
layer,  and  special  geometric  descriptions)  will  be  given. 

HYDRO  Input  Deck 


Title  Card  FORMAT  (9A10) 

TITLE  - 80  column  title. 

Foil  Component  Control  Card  FORMAT  (2F10.0) 

NC  - Number  of  foil  components  (for  HYDROSAP  this  number  is 
unity) . 

(Continued) 

*These  subroutines  are  not  required  for  the  single  element  foils  dealt  with 
in  the  present  study;  however,  they  are  retained  (and  discussed)  since 
they  constitute  a part  of  the  total  computational  program  and  may  be  of 
value  to  future  users. 


CARD  1: 


CARD  2: 
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CARD  3: 


CARD  4: 


CARD  5: 


CARD  6: 


NCPC  - number  of  computational  points  for  (each  successive) 
foil  component  (NCPC  s NCMAX  and  £ 81) 

Plot  Control  Cards  FORMAT  (4F10.0) 

(This  card  set  is  machine  and  hardware  dependent  --  a description, 
for  plotting,  is  not  given  here  since  the  routines  currently  used 
are  not  for  export.) 

Geometry  and  Panel  Control  Cards  FORMAT  (4F10.0,  4A10) 

NPP  - Number  of  pivotal  points  referenced  to  (each)  foil 
component  (0  < NPP  $ NC-1). 

SF  - Scale  actor  for  (each)  foil  component  (this  factor  con- 
verts component  geometry  and  associated  transition 
points,  as  well  as  pivotal  point  coordinates,  to  feet). 

I OP  - Panel  distribution  option  indicator  ( NCPC- 1 panels  for 
each  component). 

I0P=0.  Panels  are  distributed  so  that  regions,  on  the 
foil  surface,  with  greatest  curvature  have 
smaller  sized  panels. 

I 0P= 1 . Panels  are  distributed  by  means  of  a modified 
cosine/sinh  method. 

I0P=2.  Panels  are  distributed  in  accordance  with  (and 
the  same  as)  the  input  geometry. 

ISMO  - Maximum  number  of  coordinate  smoothing  iterations  allowed 
(for  no  smoothing  set  ISM0=0). 

CLABEL  - 40  column  component  label. 

Pivotal  Point  Coordinate  Description  FORMAT  (8F10.0) 

XP,ZP  - Coordinates  of  a pivotal  point,  referenced  to  (each) 
foil  component  (card  omitted  if  NPP=0). 

Foil  Coordinate  Count  Designation  FORMAT  (F10.0) 

NU  - Number  of  foil  upper  surface  coordinates 
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CARD  7: 


CARD  8: 


CARD  9: 


CARD  10 


CARD  11 


CARD  12 


Foil  Coordinates  (upper  surface)  FORMAT  (2F10.0) 

XU,ZU  - Foil  upper  surface  coordinates  in  (x,z)  pairs  --  start- 
ing at  the  leading  edge  and  ending  at  the  trailing  edge 
(input  NU  cards). 

Foil  Coordinate  Count  Designation  FORMAT  (F10.0) 

Nl_  - Number  of  foil  lower  surface  coordinates. 

Foil  Coordinates  (lower  surface)  FORMAT  (2F10.0) 

XI, ZL  - Foil  lower  surface  coordinates  in  (x,z)  pairs  --  start- 
ing at  the  leading  edge  and  ending  at  the  trailing  edge 
(input  NL  cards;  NOTE  THAT  NU+NL  s 200). 

Note:  CARDS  4 through  9 are  input  once,  per  foil,  since  NC=1  (in 
the  HYDROSAP  program). 

Foil  Component  Index  Card  FORMAT  (FI 0.0) 

IM  - Index  designation  for  the  main  (reference)  foil  component 
(card  omitted  when  NC=1). 

Foil  Component  Index  Card  FORMAT  (5F10.0) 

IC  - Index  for  foil  component. 

IPP  - Index  for  pivotal  point,  to  be  used,  in  "placing"  this 

component. 

ICR  - Index  for  the  foil  reference  component. 

IPPR  - Index  for  pivotal  point,  on  the  reference  component,  to 
be  used  in  "placing"  this  component. 

DELTA  - Rotation  angle  between  this  component's  coordinate  frame 
and  the  reference  component's  coordinate  frame  (in 
degrees;  clockwise  rotation  is  positive). 

Note:  CARO  11  is  repeated  NC-1  times  (card  omitted  when  NC=1). 

Angle  of  Attack  Count  Card  FORMAT  (F10.0) 

NA  - Number  of  angles  of  attack  to  be  input  (1  < NA  $ 10). 
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CARD  13:  Angles  of  Attack  Described  FORMAT  (8F10.0) 

ALPHA  - Angle  of  attack,  in  degrees  (NA  values  input). 


CARD  14:  Mach  Number  Count  FORMAT  (FI 0.0) 

NM  - Number  of  freestream  mach  numbers  to  be  input 
(1  < NM  5 10). 


CARD  15:  Freestream  Mach  Numbers  Described  FORMAT  (8F10.0) 

FSMACH  - Freestream  mach  number  (NM  values  input). 

CARD  16:  Added  Input  Information  FORMAT  (2F10.0) 

CREF  - Reference  (foil)  chord,  in  feet.  (This  value  is  used, 
internally,  to  nondimensional i ze  output  quantities.) 

RN  - Reynold's  number,  based  on  reference  chord  (CREF)  and 
freestream  conditions.  (TOTAL  RN  is  expressed  in 
MILLIONS. ) 

CARD  17:  Transition  Point  Description  FORMAT  (6F10.) 

LTRANU  - Upper  foil  surface  transition  option. 

LTRANU=0  implies  "free"  transition  ( XTRAN=ZTRAN=0. 0) . 
LTRANU=1  implies  "fixed"  transition. 

XTRANU  - Upper  surface  x-location  for  a fixed  transition. 

ZTRANU  - Upper  surface  z-location  for  a fixed  transition. 

LTRANL  - Lower  surface  transition  option. 

XTRANL  - Lower  surface  x-location  for  a fixed  transition. 

ZTRANL  - Lower  surface  z-location  for  a fixed  transition. 

Note:  CARD  17  is  repeated  NC  times  (NC=1  in  HYDROSAP). 

CARD  18:  Confluent  Boundary  Layer  Control  Card  FORMAT  (8F10.0) 

ICFL  - Confluent  boundary  layer  option  ( NC- 1 values  input; 
card  omitted  i f NC=1 ) . 

ICFL=0.  Do  not  compute  confluent  boundary  layer  for 
this  foil  component. 

(Continued) 
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ICFL=1.  Compute  confluent  boundary  layer  for  this  foil 
component. 


3 . 4 HYDRO  Output 

Output  from  the  HYDRO  subprogram,  described  below,  refers  to  the  printed 
information  produced  at  the  GSFC  computing  center  following  the  successful 
execution  of  a typical  HYDROSAP  run. 

The  printed  output,  received  from  the  computing  center,  consists  of 
all  information  requested  (both  as  special  print  output  plus  the  "usual" 
information  provided  to  the  user).  These  data  normally  appear  as  a three 
. tiered  compilation  consisting  of:  HYDROSAP  input  "echo"  data;  the  HYDRO 
subprogram’s  output  segment;  and,  the  NONSAP  output  segment.  It  is  the 
middle  segment  of  the  data  output  which  is  discussed  below. 

3.4.1  HYDRO  Print-out  Information 


This  output  segment  begins  immediately  following  the  message  titled, 

"END  ECHO".  This  message  denotes  the  terminus  for  the  HYDROSAP  input  stream, 
which  has  been  "echoed".  Should  an  echo  of  the  input  stream  not  be 
requested,  then  the  terminating  message  would  have  followed  ( immediately) 
the  namelist  (DRV)  echo. 

HYDRO  output  begins  with  an  identi f ication  message,  followed  immediately 
by  a listing  of  the  input  foil  coordinate  data  (see  section  3.3  for  notations 
and  a listing  of  the  input  quantities). 

The  input  coordinates  are  repeated  in  the  same  format  as  the  data  were 
entered  into  the  program  --  the  input  calls  for  these  coordinates  to  be 
expressed  in  inches. 

Next,  a tabulation  of  the  dimensionless,  redistributed  coordinates  -- 
with  arc  length  (SARC/C)  and  theta-equivalent  point  descriptions  (THETA) 
noted  in  parallel  vector  arrays.  The  manner  in  which  this  "redistribution" 


is  made  depends  on  the  (user  supplied)  value  given  to  variable  (IOP);  see 
section  3.3  for  discussion. 

The  SARC/C  tabulation  is  included  in  the  output  to  denote  the  surface- 
arc  length,  to  each  point  station  on  the  foil's  surface,  measured  from  the 
trailing  - along  the  lower  surface  - and  over  the  upper  surface  - (back)  to 
the  trailing  edge.  These  data  (and  the  theta-distribution  table)  are  of 
interest  to  the  aero/hydodynamists;  they  are  of  small  value  to  the  engineer. 

Following  the  table  is  a section,  in  the  printout,  headed  "Load 
Summary  Sheet  ...".  The  information  listed  here  (by  iteration  number)  is 
of  more  use  to  the  program's  monitor  than  to  the  engineer  seeking  informa- 
tion. Nonetheless,  for  completeness  and  clarity,  some  explanatory  comments 
should  be  made  regarding  the  information  contained  here. 

Recall  that  for  a Viscous/Potential  Flow  problem  the  program  must 
operate  in  an  iterative  mode,  to  obtain  a solution,  in  order  to  account  for 
boundary  layer  buildup  over  the  foil's  surface.  (In  the  event  that  a "clas- 
sical" potential  flow  solution  would  be  called  for,  there  would  be  no  need 
to  "iterate"  and  a solution  would  be  forthcoming  immediately.) 

The  tabulated  data,  presented  for  each  iteration,  lists  the  component 
number  (only  unity  will  appear  for  the  HYDROSAP  program),  followed  by  the 
"derived"  information.  These  latter  data  denote  flow-field  conditions  and 
the  force  and  moment  coefficients  which  have  been  determined  from  the  inte- 
grated pressure  - shear  distribution(s) . A header  line,  noting  physical 
conditions  for  the  evaluations  precedes  the  above  discussed  tabulation. 

Here  a Mach  number  is  displayed  since  the  program  requires  a small  input 
value,  for  this  parameter,  in  order  to  run  properly.  The  Reynold's  num- 
ber show  there  is  a value  (calculated)  representative  of  the  appropriate 
conditions  for  the  current  problem. 

The  reader  will  note  that  for  a (general  type)  viscous  flow  solution 
there  will  be  several  iterations  required  in  obtaining  a solution;  and, 
that  each  such  iterate  will  have  a tabulation  (as  described  above).  Once 


an  acceptable  solution  has  been  achieved,  the  program  does  a listing  of  the 
"final"  values  developed  within  the  various  program  segments. 

First,  there  is  a listing  headed  "Boundary  Layer  Data"  --  this  tabu- 
lates (by  surface  coordinates)  the  layer  information.  Here  a description 
of  the  laminar-turbulent  regions,  over  the  upper  and  lower  (foil)  surfaces 
is  given  (by  column  vectors  correspoinding  to  position).  Aside  from  these 
loci,  the  most  useful  information  presented  in  the  tables  will  be  the 
columns  headed  "CP"  and  "CF".  These  describe  the  local  values  for  pres- 
sure coefficient  (CP)  and  skin  function  coefficient  (CF).  The  other  data 
listed  pertain  to  a description  of  the  boundary  layer  solution  (DSTARC, 
THETMC,  H and  DELTA/C);  also  listed  are  local  flow  conditions  (ML),  and 
arc  length  (S/C)  --  as  mentioned  above. 

A notation,  following  each  tabulation,  denoting  separation,  refers  to 
the  laminar  separation;  this  and  the  reattachment  message  are  indicators 
of  other  computed  flow  field  conditions  which  exist  because  of  the  boundary 
layer  situation.  (The  reader  is  advised  (again)  that,  at  present,  we  are  un 
able  to  adequately  predict  "flow  separation"  from  the  foil's  surface.) 

Immediately  following  the  tabulation  for  boundary  layer  information  is 
the  final  load  summary  table.  Here  the  last  predicated  values  (describing 
hydrodynamic  characteristics,  stagnation  and  transition  loci,  etc.)  for  the 
foil  section  are  listed  for  the  readers'  perusal  and  subsequent  use.  The 
last  line  of  the  summary  makes  note  of  those  primary  parameters  ascertained 
from  the  analysis.  There  the  lift  coefficient  (CL),  predicted  drag  coef- 
ficient (CD),  the  Squire-Young  profile  drag  coefficient  (CDSY),  moment 
coefficients  (nose  and  quarter-chord  values)  and  the  lift-to-drag  ratio 
(CL/CD)  are  indicated. 

This  line  summary  completes  the  HYDRO  output  and  signals  a termination 
of  the  "usual"  information  displayed  from  the  subprogram. 


Following  the  HYDRO  output  section  is  a like  listing  from  the  NONSAP 
subprogram.  A discussion  of  these  output  data  will  be  found  in  section  4. 


DESCRIPTION  OF  NONSAP  - THE  NONLINEAR  STRUCTURAL  ANALYSIS  PROGRAM 

The  finite  element  program  NONSAP  (Nonlinear  Structural  Analysis  Program) 
was  originally  developed  in  1974  [11].  NONSAP  may  be  used  to  analyze  the 
static  or  dynamic  response  of  one,  two,  or  three  dimensional  structures, 
including  the  effects  of  geometric  and/or  material  nonlinearities. 

In  this  section  we  give  a brief  overview  of  NONSAP's  capabilities  and 
operation,  provide  a short  description  of  the  various  NONSAP  subroutines, 
and  describe  the  NONSAP  input  format.  However,  for  a more  complete  and 
detailed  description  of  the  NONSAP  theory  and  operation,  the  reader  is 
referred  to  the  original  program  documentation  [11,12], 

4.1  Program  Description 

The  NONSAP  program  was  developed  to  solve  static  or  dynamic,  linear  or 
nonlinear  structural  problems.  The  nonlinearities  may  be  due  to  material 
behavior  (elastic,  hyperelastic,  and  hypoelastic  material  models  are  avail- 
able) or  to  the  effects  of  large  displacements  and  large  strains.  The 
equations  of  motion  of  continuum  mechanics  governing  the  behavior  of  the 
structure  under  study  are  discretized  in  the  space  variables  through  the  use 
of  isoparametric  finite  elements.  The  equations  are  discretized  in  the 
time  variable  through  either  a Wilson-Theta  or  a Newmark-Beta  time  integra- 
tion scheme.  The  resulting  nonlinear  equations  are  solved  using  an 
incremental  method  approach,  one  which  corresponds  to  a modified  Newton 
iteration  scheme. 

The  structural  systems  to  be  analyzed  may  be  modeled  using  a variety 
of  finite  element  types;  these  are; 

• one-dimensional  truss  elements, 

• two-dimensional  plane  stress  or  plane  strain  elements, 

• two-dimensional  axi-symmetric  shell  elements,  and 
three  dimensional  solid  or  thick  shell  elements. 


Since  nonlinearities  may  be  due  to  nonlinear  material  behavior  as  well 
as  to  large  displacements  or  large  strains,  there  are  several  material  des- 
criptions available  in  the  program.  These  include: 

• an  isotropic  linear  elastic  material  model, 

• an  orthotropic  linear  elastic  material  model, 

• the  Mooney-Rivl in  model, 

• ar,  elastic-plastic  material  with  von  Mises  or  Drucker-Prager 
yield  conditions, 

• a variable  tangent  model,  and 

• a curve  description  model. 

Currently,  the  HYDROSAP  prgram  makes  use  only  of  the  one-dimensional  truss 
finite  element,  the  two-dimensional  plane  strain  element,  and  the  isotropic 
linear  elastic  material  model.  All  nonlinearities  are,  therefore,  geo- 
metric. Moreover,  HYDROSAP  uses  NO', 'SAP  routines  in  the  nonlinear  static 
solution  mode  with  the  "updated  Lagranaian"  solution  procedure.  As  a con- 
sequence, the  discussions  below  will  be  restricted  to  those  NONSAP  procedures 
and  options  which  are  currently  employed  in  HYDROSAP. 

4.1.1  Incremental  Equilibrium  Equations  for  Structural  Systems 

Using  the  well-known  and  well  documented  [11,13,14]  method  of  finite 
elements,  the  incremental  nodal  point  equilibrium  equations  for  an  assem- 
blage of  finite  elements  take  the  form 

M Li ( t+A t)  + C u ( t+At ) + K( t ) au  = R(t+At)  - F(t)  (1) 

where 


M = the  constant  mass  matrix, 

C = the  constant  damping  matrix, 

K(t)  = the  tangent  stiffness  matrix  at  time  t , 

R(t)  = the  external  applied  load  vector  at  time  t , 

F(t)  = the  nodal  point  force  vector  (of  stress-equivalent 

internal  reactions)  at  time  t , 
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u(t)  = the  vector  of  nodal  point  displacements,  at  time  t , 

au  = the  vector  of  nodal  point  displacement  increments  from 

time  t to  t+At:  u(t+At)  - u(t)  , 
t = time,  and 

At  = time  increment. 

Of  course,  in  a static  analysis  M=C=0  , while  time  is  regarded  only  as  a 
loading  parameter  used  for  applying  the  external  load  in  "quasistatic"  load 
increments.  Also,  in  a nonlinear  analysis,  the  solution  of  (1)  yields 
only  approximate  solution  increments,  au  . In  order  to  improve  the  solu- 
tion accuracy  and  to  prevent  the  buildup  of  solution  instabilities,  equilib- 
rium iterations  may  be  used  at  selected  time  steps. 

In  the  nonlinear  static  case  an  equation  system  of  the  form 

K(t)  Uu)i+i  = R(t+At)  - Fi  i=0,l ,2 , • • • (2) 

is  solved  by  iteration.  Here  Fq  = F(t)  and  F^  is  the  vector  of  stress- 
equivalent  internal  reactions  corresponding  to  the  displacement  vector 
u^t+At)  = ui _ 1 ( t-*-At ) + (au)^  . The  solution  of  (2),  by  iteration,  is 
clearly  a modified  Newton  method  for  (1). 

4.1.2  NONSAP  Solution  Process 

The  complete  NONSAP  solution  process  consists  of  three  distinct  stages: 
input,  assembly,  and  step-by-step  solution  (see  Figure  4-1). 

Input  Stage 


In  this  stage,  control  and  nodal  point  input  are  read  and/or  generated 
by  the  program;  and  equation  numbers  for  the  active  degrees  of  freedom  at 
each  nodal  point  are  established.  In  addition,  externally  applied  load 
vectors  are  established  and  then  stored  on  tape.  Finally,  the  element 
data  are  ready  and/or  generated,  element  connectivity  arrays  are  established, 
and  all  element  dataare  then  stored  on  tape. 
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Figure  4-1.  NONSAP  Flowchart 

Program  Operation  In  Nonlinear 
Static  Solution  Mode 
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Assembly  Stage 


Here  the  linear  and  effective  linear  structural  stiffness  matrices 
are  assembled  and  stored  on  tape,  as  are  the  mass  and  damping  matrices 
(for  a dynamic  analysis).  Matrices  are  stored  in  compacted  form  as  one- 
dimensional arrays.  Only  the  elements  below  the  "skyline"  of  a matrix  are 
processed,  thereby  reducing  storage,  input/output,  and  computational  costs. 

Solution  Stage 

In  a nonlinear  static  analysis  (the  solution  mode  for  which  NONSAP 
routines  are  used  by  HYDROSAP)  damping  and  mass  effects  are  neglected. 

Time  serves  as  a loading  parameter,  only,  and  the  input  time  step  deter- 
mines the  load  increment  to  be  applied  to  the  structure  during  each  solu- 
tion step.  The  linear  stiffness  matrix  which  was  previously  assembled 
and  stored  on  tape,  is  updated  at  each  solution  step.  This  is  accomodated 
by  updating  the  stiffness  matrices  of  the  nonlinear  elements  to  form  a cur- 
rent tangent  stiffness  matrix.  The  load  step  interval  at  which  this  up- 
date occurs  is  one  of  the  NONSAP  control  inputs. 

Moreover,  within  each  load  step,  solution  accuracy  may  be  improved  by 
equilibrium  iterations.  The  interval  of  load  steps  at  which  an  equilib- 
rium iteration  is  to  be  allowed  is  also  a NONSAP  control  input. 

A solution  to  the  incremental  linear  equations  at  each  load  step  is 
performed  by  the  linear  equation  solver,  C01S0L;  a routine  which  performs 
Gauss  eliminations  on  the  positive  definite  symmetrical  system  of  equations. 
This  algorithm  takes  advantage  of  the  banded  and  sparse  structure  of  the 
stiffness  equations  by  not  processing  elements  outside  a matrix  skyline, 
since  they  remain  zero  throughout  the  computation.  The  algorithm  con- 
sists of  an  LDL*  decomposition  of  the  tangent  stiffness  matrix  into  lower 
triangular,  L , and  diagonal,  D , matrices,  followed  by  a reduction  and 
back  substituion  of  the  load  vector. 


4.1.3  NONSAP  Element  Library 


In  this  section  we  will  briefly  describe  only  those  NONSAP  element 
types  which  are  used  by  HYDROSAP. 

The  One-Dimensional  Truss  Element 


A one-dimensional  truss  element,  which 
may  be  located  in  a three-dimensional  space, 
is  used  in  HYDROSAP  to  model  the  thin  skin 
which  is  on  the  exterior  of  some  cable  fair- 
ings. This  element  has  two  nodes;  it  is 
assumed  to  have  a constant  cross-sectional 
area  (which  corresponds  to  the  skin  thickness 
in  the  HYDROSAP  plane  strain  analysis)  and  to 
undergo  small  strains  and  large  displacements. 


Z 


Figure  4-2.  A One-Dimensional 
Truss  Element 


The  Two-Dimensional  Plane  Stress/Plane  Strain  Element 


An  isoparametric  plane  strain  element, 
having  a variable  number  of  nodes,  is  used  to 
model  the  interior  (cross-section)  of  the  cable 
fairing.  This  element  may  have  from  three  to 
eight  nodes;  any  of  the  nodes  numbered  from  five 
through  eight  on  the  sketch  may  be  omitted.  A 
three  node  element  is  obtained  by  setting  the  co- 
ordinates of  two  nodes,  for  a four  node  quadri- 
lateral, equal  to  each  other.  In  a plane  strain 
analysis  the  out-of-plane  element  thickness  is  as- 
sumed to  be  unity. 

4.2  Subroutine  Descriotions 


Figure  4-3.  A Two-Dimensional 
Plane  Strain  Isoparametric 
Quadri lateral  Element 


In  this  section,  a brief  description  of  the  various  NONSAP  subroutines 
is  provided.  In  addition,  a structure  map  of  the  NONSAP  module  can  be 
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found  in  Appendix  A. 


ADDBAN  Assembles  upper  triangular  element  stiffness  into  a compact  global 
stiffness  matrix. 

ADDCM*  Adds  nodal  masses  and  dampers  to  the  global  mass  matrix  or  to 
the  global  effective  stiffness  matrix. 

ADDDI  Adds  displacement  increments  to  the  working  displacement  array. 

ADDMA*  Adds  element  masses  to  the  global  lumped  mass  array. 

ADDRES  Calculates  addresses  for  diagonal  elements  in  a banded  matrix 

whose  column  heights  are  known. 

ASSEM  Assembles  matrices  needed  in  a solution  and  dumps  them  to  mass 
storage. 

BANDET*  Performs  triangular  factorization  and  determinant  calculations  for 
a banded  matrix. 

CACCHY  Converts  Piola-Kirchoff  stresses  to  Cauchy  stresses. 

CDMOD*  Provides  the  curve  description  non-linear  material  model  (for 
two-dimensional  elements). 

CM0D3D*  Provides  the  curve  description  non-linear  material  model  (for 
three-dimensional  elements). 

COLHT  Computes  column  heights  for  a banded  sparse  matrix. 


‘These  subroutines  are  not  required  for  the  single 
in  the  present  study;  however,  they  are  retained 
they  constitute  a part  of  the  total  computational 
value  to  future  users. 


element  foils  dealt  with 
(and  discussed)  since 
program  and  may  be  of 


47 


COLSOL  The  simultaneous  linear  equation  in-core  solver  routine.  It 
uses  compacted  storage  and  a column  reduction  scheme. 


DCRACK*  This  routine  modifies  the  stress-strain  model  to  account  for 
material  cracking. 

DERIQ  A routine  which  evaluates  the  strain-displacement  matrix  for 
quadrilateral  elements  with  axisymmetric  geometry. 


DERIQ3  Evaluates  the  strain-displacement  matrix  for  the  8 to  21  node 
curvilinear  hexahedron. 


DRUCK*  Provides  element  stresses  for  the  Drucker-Prager  non-linear 
material  model . 

ELCAL  Makes  the  calls  to  appropriate  element  routines  for  reading, 

limited  checking,  generating,  and  storing  of  finite  element  data. 

ELEMNT  Provides  for  calls  to  the  appropriate  one-,  two-,  or  three- 
dimensional  element  initialization  routines. 


ELPAL*  An  algorithm  to  generate  the  elasto-plastic  non-linear  material 
model . 


ELT2D3*  The  2-D  material  model  #3:  Variable  Tangent  Moduli  Model 


ELT2D4*  The  2-D  material  model  #4  or  #5:  Curve  Description  Model 


ELT2D6*  The  2-D  material  model  #6:  Von  Mises  Elastoplastic  Model. 


ELT2D7*  The  2-D  material  model  t> 7:  Drucker-Prager  Elastoplastic  Model 


*These  subroutines  are  not  required  for  the  single  element  foils  dealt  with 
in  the  present  study;  however,  they  are  retained  (and  discussed)  since 
they  constitute  a part  of  the  total  computational  program  and  may  be  of 
value  to  future  users. 


r4 
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ELT2D8*  The  2-D  material  Node!  #8:  Incompressible  Elastic  Model 
ELT2D9*  A dummy  routine  for  any  user-provided  non-linear  material  model. 


ELT3D3*  The  3-D  material  model  #3:  Curve  Description  Non-linear  Model. 


ELT3D4* 
ELT3D5* 
ELT3D6* 
EL2D1 0* 
EL2D11* 
EL2D12* 


Dummy  material  models  (user-supplied). 


EQUIT  An  iteration  scheme  used  to  achieve  dynamic  equilibrium. 


ERROR  This  routine  prints  error  message  if  storage  requirements  exceed 
program  capaci ty . 

FREQS*  This  routine  finds  lowest  frequencies  and  mode  shapes  for  the 
linearized  structure. 

FUNCT*  A subprogram  which  computes  the 

• interpolation  functions  and  their  derivatives  for  the 
curvilinear  isoparametric  hexahedron  (8  to  21  nodes), or 
the 

• Jacobian  of  element  mapping  function. 


FUNCT2  A routine  to  compute  the 

• interpolation  functions  and  their  derivatives  for  the  curvi- 
linear isoparametric  quadri lateral  (4  to  8 nodes),  or  the 

• Jacobian  of  element  mapping  function. 


♦These  subroutines  are  not  required  for  the  single  element  foils  dealt  with 
in  the  present  study;  however,  they  are  retained  (and  discussed)  since 
they  constitute  a part  of  the  total  computational  program  and  may  be  of 
value  to  future  users. 
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ICDMOD*  An  initializing  routine  for  CDMOD. 


ICM0D3*  An  initializing  routine  for  CM0D3D. 


IDRUCK*  An  initializing  routine  for  DRUCK. 


IELPAL*  An  initializing  routine  for  ELPAL. 


INITAL  Initializes  the  displacement,  velocity,  and  acceleration  arrays. 

INI TWA*  Initializes  the  working  array  used  with  the  2-D  non-linear  material 

models. 


INPUT  A routine  which 


reads  and  prints  nodal  point  data, 
generates  missing  nodal  data,  and 

determines  equation  numbers  associated  with  each  node. 

INTWA3*  Initializes  the  working  array  used  with  the  3-D  non-linear  material 
models . 

IVTMOD*  An  initializing  routine  for  VTMOD. 

LOADEF  An  algorithm  to  calculate  effective  applied  loads  (excluding  non- 
linear contributions). 

LOADER  (A  BTS-supplied  HYDROSAP  routine)  converts  a pressure  distribution 
to  consistent  nodal  loads  and  assigns  appropriate  load-profile 
multipliers  for  each  direction  at  every  node. 


‘'These  subroutines  are  not  required  for  the  single  element  foils  dealt  with 
in  the  present  study;  however,  they  are  retained  (and  discussed)  since 
they  constitute  a part  of  the  total  computational  program  and  may  be  of 
value  to  future  users. 
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LOADS 


This  subroutine 


• reads  load  profiles  and  interpolates  functional  values  at 
requested  points  in  time, 

• reads  nodal  load  factors,  and 

• calculates  noaal  loads  for  each  time  step  and  stores  the 
results  on  Tape  3. 

MATBAR  Computes  stresses  for  the  1-D  truss  elements. 


MATRT2  Prints  material  properties  for  the  2-D  elements. 

MATVJRT*  Prints  material  properties  for  the  3-D  elements. 

MAXAEX*  Shifts  elements  of  a working  array  during  the  eigenvalue  analysis. 

MAXMIN  An  algorithm  to  calculate  principal  stresses. 

MIDEP*  A routine  which  forms  the  elasto-plastic  material  matrix. 


MLTPLY*  An  algorithm  to  compute  A=B*C  , where  B is  a matrix  (stored 

in  compacted  form)  and  A and  C are  vectors. 

MOONEY*  An  algorithm  to  describe  an  incompressible  nonlinear  elastic 

material  (Mooney-Ri vl in  description  in  a state  of  plane  stress). 

MULT  An  algorithm  to  compute  A=B*C  , where  B is  a matrix  (stored 

in  compacted  form)  and  A and  C are  vectors. 

UDLLDS  (A  BTS-supplied  HYDROSAP  routine.)  Determines  consistent  nodal 

loads  from  the  fairing  surface  geometry  and  its  pressure  distribution. 


*These  subroutines  are  not  required  for  the  single  element  foils  dealt  with 
in  the  present  study;  however,  they  are  retained  (and  discussed)  since 
they  constitute  a part  of  the  total  computational  program  and  may  be  of 
value  to  future  users. 
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NEWDAV  This  routine  calculates  the  new  displacements,  velocities,  and 
accelerations  at  time  t+At. 

NODMAS*  A routine  that 

• reads  concentrated  nodal  masses  and  dampers,  and 

• computes  nodal  mass  and/or  damping  vectors  and  stores  them 
on  Tape  1 1 . 

NONSAP  The  main  program  for  the  NONSAP  module,  in  HYDROSAP.  This  routine 
controls  NONSAP  execution. 

NONSPI  A subprogram  to  read  MONSAP  input  files. 

OPCOEF  This  subprogram  calculates  the  coefficients  for  the  time-integra- 
tion operators  used  in  the  Wilson  Theta  Method  or  the  Hewmark 
Beta  Method. 

PNORM*  An  algorithm  to  compute  a matrix  pseudo-norm. 

POSINV  A subprogram  used  to  invert  a positive  definite  matrix. 

PRAGER*  Forms  the  el  as to-pl astic  material  matrix  for  the  Drucker-Prager 
material  model . 

QUADM  Calculates  the  consistent  mass  matrix  for  a quadrilateral  element. 

QUADM3*  Calculates  the  consistent  mass  matrix  for  a curvilinear  hexahedron. 

QUADS  Determines  the  element  stiffness  matrix  for  quadrilateral  axisym- 

metric  or  plane-stress/plane-strain  2-D  elements. 


’•'These  subroutines  are  not  required  for  the  single 
in  the  present  study;  however,  they  are  retained 
they  constitute  a part  of  the  total  computational 
value  to  future  users. 


element  foils  dealt  with 
(and  discussed)  since 
program  and  may  be  of 
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QUADS 3* 


Determines  the  element  stiffness  matrix  for  3-D  hexahedral  curvi- 
1 inear  elements . 


READID  Reads  the  ID  array  from  Tape  8. 

RUSS  Reads  and  generates  1-D  truss  element  data;  also  assembles  the 

1- D  truss  elements. 

SECANT*  This  routine  computes  eigenvalues/eigenvectors  for  a modal 
analysis. 

ST  I ME  An  algorithm  which  computes  CPU  time  (seconds)  remaining  in  a job 
time  estimate. 

STRESS  Calls  subroutine  ELEMNT  for  stress  computations. 

STSTL  Generates  the  global  stress-strain  law  for  isotropic/orthotropic 

2- D  elements. 

STSTN*  Generates  the  stress-strain  law  and  stresses  for  nonlinear  material 
models . 

STST3L*  Generates  the  global  stress-strain  law  for  isotropic/orthotropic 

3- D  elements. 

STST3N*  Generates  stresses  for  all  3-D  models,  and  the  stress-strain  law 
for  nonlinear  material  models. 

TDFE  This  routine  controls  the  generation  and  assembly  of  2-D  elements. 

THDFE*  This  routine  controls  the  generation  and  assembly  of  3-D  elements. 


*These  subroutines  are  not  required  for  the  single  element  foils  dealt  with 
in  the  present  study;  however,  they  are  retained  (and  discussed)  since 
they  constitute  a part  of  the  total  computational  program  amd  may  be  of 
value  to  future  users. 
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THREDM* 


This  subroutine 


• allocates  storage  for  the  3-D  elements,  and 
calls  THDFE. 

TODMFE  This  routine 

• allocates  storage  for  the  2-D  elements,  and 
calls  TDFE. 

TTIME  An  algorithm  which  computes  elapsed  CPU  time  (in  seconds)  from 
last  call  to  STIME. 

TRUSS  This  routine 

allocates  storage  for  the  1-D  elements,  and 
calls  RUSS. 

VTMOD'  The  variable  tangent  moduli  material  model. 

WRITE  A subroutine  to  print  displacements,  velocities  and  accelerations. 

WR1FRQ  A subroutine  to  print  frequencies  and  mode  shapes  in  an  eigen- 

system  analysis. 

4.3  NQNSAP  Input 

This  section  details  the  contents  of  the  NONSAP  input  file.  With 
KYDROSAP  operating  in  the  SHRTIN  = TRUE  mode,  this  file  is  built  auto- 
matically. However,  for  specialized  applications  the  user  may  wish  to 
construct  this  file  himself.  We  will  omit  discussing  the  NONSAP  control 


♦These  subroutines  are  not  required  for  the  single  element  foils  dealt  with 
in  the  present  study;  however,  they  are  retained  (and  discussed)  since 
they  constitute  a part  of  the  total  computational  program  and  may  be  of 
value  to  futureusers. 
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MMM 


card  groups  (such  as  those  governing  three-dimensional  or  dynamic  analysis) 
which  are  not  applicable  to  a nonlinear,  static,  plane  strain  analysis. 


A complete 
where  [12]. 

CARD  1: 

CARDS  2: 

Card  2.1 


discussion  of  the  NONSAP  input  format  is  available  else- 


NONSAP  Input  Deck 


Heading  Card  FORMAT  (10A8) 

HED  - 80  column  title. 

Master  Control  Cards  FORMAT  (15,  611,  14, 

315,  2F10.0,  15) 


NUMP 

- 

Total  number  of  nodal  points; 

= 0,  program  stops. 

I DOF ( 1 ) 

Master  X-translation  code 

= 0,  admissible  degree  of  freedom 
= 1,  deleted  degree  of  freedom. 

I DOF ( 2 ) 

- 

Master  Y-translation  code. 

1D0F(3) 

- 

Master  Z-translation  code. 

I DOF ( 4 ) 

- 

Master  X-rotation  code. 

I DO  F ( 5 ) 

- 

Master  Y-rotation  code. 

I DOF ( 6) 

- 

Master  Z-rotation  code. 

NEGL 

- 

Number  of  linear  element  groups; 

= 0,  all  elements  are  nonlinear. 

NEGNL 

- 

Number  of  nonlinear  element  groups 

= 0,  all  elements  are  linear. 

MODEX 

Solution  mode; 

= 0,  data  check  only 

= 1 , execution 

= 2,  restart. 

NSTE 

- 

Number  of  solution  time  steps. 

DT 

- 

Time  step  increment. 

TSTART 

- 

Time  at  solution  start. 

IPRI 

- 

Output  print  interval. 
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Card  2.2 


Card  2.3 


Card  2.4 


Card  2.5 


F0RI1AT  (415) 

IMASS  - Static  or  dynamic  analysis  control  flag; 

= 0,  static  analysis  (use  only  this  value  in 
a full  HYDROSAP  analysis) 

= 1,  dynamic  analysis  with  lumped  (diagonal) 
mass  matrix 

= 2,  dynamic  analysis  with  consistent  mass 
matrix. 

I DAMP  - Damping  control  flag; 

= 0,  no  damping  (use  only  this  value  in  a full 
HYDROSAP  analysis) 

= 1,  Rayleigh  damping  assumed. 

IMASSN  - Number  of  concentrated  nodal  masses. 

IDAMPN  - Number  of  concentrated  nodal  dampers. 

FORMAT  (15) 

IEIG  - Eioensystem  control  flag; 

= 0,  no  eigensystem  solution 
= 1,  find  lowest  frequencies  and  mode  shapes. 

FORMAT  (415,  E10.4) 

1SREF  - Number  of  time  steps  between  effective  stiff- 

ness matrix  updates. 

NUMREF  - Number  of  allowable  stiffnes  matrix  updates 
within  each  time  step. 

IEQUIT  - Number  of  time  steps  between  equilibrium 
i terations . 

ITEMAX  - Maximum  number  of  equilibrium  iterations  permit- 
ted . 

RTOL  - Relative  tolerance  used  to  measure  equilibrium 

convergence. 

FORMAT  (110,  2 FI  0.0) 

IOPE  - Time  integration  method; 

= 1,  Wil son-Theta  method 
= 2,  Newmark-Beta  method. 
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(Continued ) 


Card 


Card 


CARDS  3: 


OPVAR(I)  - First  integration  parameter. 

0PVAR(2)  - Second  integration  parameter. 


Note:  This  card  is  read  but  ignored  in  a static  analysis 
( IMASS=0) . 


2.6 

NPB 


1DC 


1VC 
I AC 


FORMAT  (415) 

- Number  of  blocks  for  displacement/velocity/ 
acceleration  output; 

= 0,  print  all  nodal  components. 

- Displacement  print  code 
= 0,  no  print 

= 1 , print. 

- Velocity  print  code. 

- Acceleration  print  code. 


2.7 


IPNODE ( 1 .1  )- 
I PNODE (2,1  )- 
I PNODE ( 1 ,2)- 


FORfiAT  (1615) 

First  node  of  printout  block  no.  1. 
Last  node  of  printout  block  no.  1. 
First  node  of  printout  block  no.  2. 


Note:  Leave  this  card  blank  if  NPB=0. 

Nodal  Point  Data  FORMAT  (A1 , 14,  A1 , 14, 

515,  3F10.0,  15) 

CT  - Coordinate  system  for  this  node; 

= (blank),  cartesian  (X,Y,Z) 

= X,  X-cylindrical . 

N - Node  number  of  this  node  (1<N<NUMNP) 
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PSF 


ID(1,N) 

ID(2,N) 

ID(3,N) 

ID(4,N) 

I D ( 5 , N ) 

I D ( 6 , N ) 

X(N) 

Y(N) 

Z(N) 

KN 


- Print  suppression  flag; 

= (blank),  no  suppression 

= A,  suppress  ordered  list  of  node  coordinates 
= B,  suppress  list  of  equation  numbers 
= C,  both  A and  B above. 

- X-translation  boundary  code. 

- Y-translation  boundary  code. 

- Z-translation  boundary  code. 

- X-rotation  boundary  code. 

- Y-rotation  boundary  code. 

- Z-rotation  boundary  code. 

- X (or  Z)  coordinate  (inches). 

- Y (or  R)  coordinate  (inches) 

- Z (or  9)  coordinate  (inches  or  radians) 

- Node  number  increment  for  automatic  nodal  data 
generation 


Note:  Input  one  card  per  node  being  defined. 


CARDS  4: 

Card  4.1 


Applied  Loads  Data  FORMAT  (315) 

NLOAD  - Number  of  nodal  load  cards. 

NLCUR  - Number  of  load  profile  time  functions. 

NPTM  - Maximum  number  of  points  used  to  describe 

any  one  load  curve 


Cards  4.2  Load  Curve  Data 


Card  4.2.1  Control  Data  FORMAT  (215) 

NTF  - Time  function  number  (1 sNTFjNLCUR) 

ilPTS  - Number  of  ordered  pairs  defining  the  time 

function  (2fNPTSfNPTM) . 

Card  4.2.2  (T,F(T) ) Pairs  FORMAT  (8F10.0) 

T I MV ( 1 ) - Time,  T , at  point  1 . 

RV(1)  - Function  value,  F(T)  , at  point  1. 

T I MV ( 2 ) 
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(Continued) 


Card  4.2 


CARD  5: 


CARDS  6: 


RV(2) 

TIMV(NPTS) 

RV(NPTS) 

.3  Nodal  Loads  Data  FORMAT  (315,  F10.0) 

Note:  Skip  this  card  if  NLOAD=0  , otherwise  input 
NLOAD  cards. 

NOD  - Node  number  to  which  this  load  is  applied. 

ID1RN  - Degree  of  freedom  (direction)  in  which  this 
load  is  applied; 

= 0,  solution  terminated 
= 1,  X-translation 
= 2,  Y-translation 
= 3,  Z-translation. 

FAC  - Function  multiplicative  factor  used  to  scale 

F(T)  for  this  load. 

Rayleigh  Damping  Specification  FORMAT  (2F10.0) 

(Omit  this  card  if  I DAMP  = 0.) 

A DAMP  - Rayleigh  damping  coefficient  a . 

BDAMP  - Rayleigh  damping  coefficient  8 . 

Concentrated  Nodal  Masses  FORMAT  (110,  6F10.0) 

(Omit  these  cards  if  IMASSN  = 0 , otherwise  input 
IMASSN  cards . ) 

N - Node  number. 

XMASS(l)  - X-direction  mass. 

XMASS(2)  - Y-direction  mass. 

XMASS(3)  - Z-direction  mass. 

XMASS(4)  - X-rotational  mass. 

XMASS(5)  - Y-rotational  mass. 

XMASS(6)  - Z-rotational  mass. 


59 


CARDS  7: 


CARDS  8: 

Card  8.1 


Cards  8. 


CARDS  9: 


Concentrated  Nodal  Dampers  FORMAT  (110,  6F10.0) 

(Omit  these  cards  if  IDAMPN  = 0 , otherwise  input  IDAMPN 
cards . ) 

N - Node  number. 

XDAMP(l)  - X-di recti  on  damper. 

XDAMP(2)  - Y-direction  damper. 

XDAMP(3)  - Z-direction  damper. 

XDAMP(4)  - X-rotation  damper. 

XDAMP(5)  - Y-rotation  damper. 

XDAMP(6)  - Z-rotation  damper. 


Initial  Conditions 


FORMAT  (15) 


ICON  - Initial  condition  flag 

= 0,  zero  initial  conditions  generated 
automatically 

= 1,  non-zero  initial  conditions  read  from 
input. 


DISP(l) 


DISP(NEQ) 

VEL(l) 


FORMAT  (6E1 2. ) 

Initial  displacement  vector,  NEQ  = number 
of  equations. 

J 

1 Initial  velocity  vector  (omit  if  IMASS  = 


VEL(NEQ)  J 
ACC ( 1 ) 

Initial  acceleration  vector  (omit  if  IMASS  = 
•*  0). 

ACC (NEQ) 


Truss  Elements 

(omit  group  if  no  truss  elements  present.) 
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(Continued) 


Cards  9. 


Cards  9. 


Card  9.2 


Card  9.2 


Card  9.3 


Element  Control  Card  FORMAT  (2014) 

NPAR(l)  - = 1 

NPAR ( 2 ) - Number  of  truss  elements  in  this  group  {>,  1) 

NPAR(3)  - Type  of  analysis: 

= 1,  material  nonlinearity  only 
= 2,  updated  Lagrangian 

NPAR(4)  ) 


- not  used 


NPAR(14)  J 

NPAR(15)  - Material  model  number; 

= 1 , linear  elastic 
= 2,  nonlinear  elastic 

NPAR (16)  - Number  of  different  sets  of  material  section 

properties  (5  1 ) 

NPAR(17)  - Number  of  material  model  constants  per  set; 

= 0,  if  NPAR ( 1 5)  = 1 
5 4,  if  NPAR( 15)  = 2 

Linear  Elastic  Material/Section  Property  Cards 
(Omit  this  group  is  NPAR(15)  / 1,  otherwise  read 
NPAR(16)  sets  of  cards.) 

. 1 Material  Number  FORMAT  (15) 

N - Material /Section  number  (1  <N<I!PAR(1  C) ) 

FORMAT  (4F10.0) 

- Young' s modul us  (psi ) 

2 

- Cross-section  area  (in  ) 

- Mass  density  (slugs) 

- Initial  axial  strain  (in/in). 

Nonlinear  Elastic  Material /Section  Property  Cards 
(Omi tted . . .not  applicable  for  the  linear  material  models 
assumed  in  HYDROSAP.) 


.2  Properties 
E(N) 
AREA(N) 

DEN (N) 
STRAl(N) 
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CARDS  10: 

Card  10.1 


2/D  Continuum  Elements 

Element  Group  Control  Card  FORMAT  (2014) 

NPAR(l)  - = 2 

NPAR(2)  - Number  of  elements  in  this  group  (>  1) 
NPAR(3)  - Type  of  analysis; 

= 1,  material  nonlinearity  only 
= 2,  total  Lagrangian 
= 3,  updated  Lagrangian 


NPAR(4 ) 
NPAR ( 5 ) 


NPAR(6) 
NPAR ( 7 ) 

NPAR(8) 
NPAR ( 9 ) 

NPAR(IO) 

NPAR(ll) 
NPAR(1 2) 

NPAR(13) 

NPAR( 1 4 ) 
NPAR ( 1 5) 


- Not  used 

- Element  type 

= 0,  axisymmetric 
= 1 , plane  strain 
= 2,  plane  stress 

- Not  used 

- Maximum  number  of  nodes  in  any  element 
(4<NPAR(7)S8) 

| Not  used 

- Numerical  integration  order  to  be  used  in 
Gauss  quadrature  (2<;NPAR(10)<4) 


Not  used 


- Number  of  stress  output  location  tables; 

= 0,  print  stresses  at  integration  points 

- Not  used 

- Material  model  number  (1 sNPAR(l 5)sl2) 

= 1,1  inear  isotropic 

= 2,  linear  orthotropic 
= 3,  variable  tangent  moduli 
= 4,  curve  description  model 
= 5,  curve  description  with  tension  cutoff 
= 6,  elastic-plastic  (von  Mises) 

= 7,  elastic-plastic  (Drucker-Prager ) 
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(Continued) 


NPAR( 16) 
NPAR ( 1 7 ) 


NPAR( 18) 


= 8,  incompressible  nonlinear  elastic 
(Mooney-Ri vl in) 

= 9,  user-supplied  material  model 

- Number  of  sets  of  material  properties 

- Number  of  material  constants  per  property  set: 

= 0,  if  NPAR( 1 5)<9 

= 1,  if  NPAR ( 1 5)>9 

- Dimension  of  storage  array  required  for  element 
history; 

= 0,  if  NPAR( 1 5)<9. 


Cards  10.2  Material  Pi operty  Data 

Input  NPAR ( 16)  sets  of  cards 


Card  10.2.1  Material  Number 


- Material  number 


DEN(N) 


FORMAT  (15,  F10.0) 


Mass  density  of  material  (slugs) 


Card  10.2.2  Material  Property  Cards  FORMAT  (2F10.0) 

(The  format  is  given  for  NPAR ( 1 5 ) = 1 only.  Other  material 
models  are  not  used  in  HYDROSAP.) 

PROP ( 1 , N ) - Young's  modulus  (psi) 

PR0P(2,N)  - Poisson's  ratio 

Card  10.3  Stress  Output  Tables 

(Not  applicable. . .omit. ) 

Card  10.4  Element  Data  Cards  FORMAT  (15,  13,  12, 

2F10.0,  1015) 

M - Element  number  (within  this  group) 

IEL  - Number  of  nodes  used  to  describe  this  element 

IPS  - Number  of  the  stress  table  to  be  used  for  this 

element; 

« 0,  no  stress  output  for  this  element 
BET  - Material  angle,  b , used  with  the  linear 

orthotropic  material  model 

6 3 

(Continued) 


THIC  - Element  thickness  (applies  to  plane  stress 

only)  (inches). 

MTYP  - Material  property  set  number  assigned  to  this 

el ement 

KG  - Node  generation  parameter  for  missing  elements 

N0D(1)  - Global  node  number  of  element  nodal  point  1 

N0D(8)  - Global  node  number  of  element  nodal  point  8. 

CARDS  11:  3/D  Continuum  Elements 

(Not  used  in  HYDROSAP. . .omitted. ) 

CARDS  12:  Eigensystem  Solution 

(Not  used  in  HYDROSAP. . .omitted. ) 

A . 4 NONSAP  Output 

In  this  section  we  discuss  the  output  format  of  the  structural  analysis 
module.  The  output  for  a sample  analysis  is  discussed  in  detail  in  section 
5.  The  NONSAP  output  (written  to  an  output  unit  defined  by  HYDROSAP  vari- 
able IWRT0)  may  be  divided  into  two  broad  classes:  problem  definition  out- 
put and  displacement  output. 

4.4.1  Problem  Definition  Output 

Problem  definition  output  follows  the  problem  title  and  consists  of  a 
summary  of  user-supplied  --  as  well  as  program-generated  --  control  informa- 
tion for  the  run;  nodal  point  data;  equation  number  assignments;  nodal 
loads  data;  and  element  definition  data.  The  list  of  control  information 
for  the  run  is  for  the  most  part  self-explanatory.  The  nodal  point  data 
gives  the  boundary  condition  codes  ( 1 constrained , 0=free)  for  each  nodal 
degree  of  freedom  and  the  cartesian  coordinates  for  each  node.  The  equa- 
tion number  data  identifies  the  equation  associated  with  each  active  degree 
of  freedom  in  the  structure. 
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The  nodal  loads  are  defined  by  a combination  of  (non-dimensional)  load 
profile  time  function(s)  and  specific  load  curve  multipliers  which  correspond 
to  each  structural  degree  of  freedom  (node  and  direction).  When  NONSAP  is 
being  used  as  part  of  HYDROSAP,  in  the  "full"  iteration  mode,  the  load 
curve  multipliers  are  determined  automatically  from  the  fairing's  (deformed) 
geometry  and  current  surface  pressure  distribution.  The  load  curve  multi- 
pliers represent  consistent  finite  element  concentrated  loads  which  are 
work-equivalent  to  the  fairing's  surface  pressure  distribution.  The  load 
profile  time  function  defined  by  the  HYDROSAP  short  form  input  pre-processor 
is  shown  in  Figure  4-4.  Note  that  the  convention  used  in  HYDROSAP  is  to 
relate  the  value  of  the  NONSAP  variable  TIME  to  the  HYDROSAP  "outer"  itera- 
tion number:  at  the  start  of  a NONSAP  computation,  in  an  outer  iteration 
ITER, 


TIME  - FL0AT( ITER-1 ) 

Thus,  for  example,  during  HYDROSAP  outer  iteration  number  five,  TIME 
increases  from  TIME=4.0  to  TIME=5.0  . In  the  current  HYDROSAF  version 
the  time  is  incremented  to  steps  of  size  DT=1.0  , except  during  outer 
iteration  I TER= 1 . During  this  first  outer  iteration  TIME  is  incremented 
in  NSTE  (sHYDROSAP  control  variable  NSTEPS)  steps  of  size  DT=1/NSTE  . 

Due  to  the  shape  of  the  load  curve  profile  (Figure  4-4),  this  procedure  has 
the  effect  of  loading  the  original  fairing  with  the  first  HYDRO  pressure 
distribution  in  "NSTE"  quasistatic  load  steps. 


r 
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The  element  definition  data  is  reported  in  element  groups.  Elements 
within  a group  are  of  the  same  type  (1-D  Truss,  2-D  quadrilateral,  3-D 
Hexahedral),  the  same  degree  of  nonlinearity,  the  same  sub-type  (axisym- 
metric,  plane  strain,  plane  stress);  and,  have  the  same  type  of  material 
model.  There  may  be  several  material  property  sets  for  each  group,  how- 
ever. For  example,  fairing  CORE  and  TAIL  elements  are  2-D  quadrilateral, 
plane  strain,  updated  Lagrangian,  isotropic  elements  which  have  differing 
elastic  properties.  Elements  are  numbered  within  each  group,  and  a 
"connect! vi ty"  table  is  printed  which  identifies  the  global  nodes  defin- 
ing each  element. 

<5.4.2  Displacement  Output 

Following  the  problem  definition  data,  the  program  prints  the  initial 
displacements  for  each  node,  and  then  the  computed  nodal  displacements 
after  every  IPRI  (z  HYDROSAP  control  variable  I PR  I NT ) steps  during  the 
solution  process.  Note  that  for  HYDROSAP  calculations  all  displacements 
are  zero  except  those  in  the  Y-Z  plane.  This  is  the  reference  plane  for 
a NONSAP  plane  strain  analysis.  The  Y-direction  is  defined  by  the  unde- 
formed fairing  chord,  and  the  Z-direction  is  transverse  to  the  chord. 

(Note  that  the  corresponding  HYDRO  coordinate  system  is  described  in  the 
X-Z  plane. ) 
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5.  SAMPLE  ANALYSIS 


To  illustrate  the  capabilities  of  the  HYDROSAP  software  system, 
the  results  of  a sample  analysis  are  discussed  in  this  section.  A 
hypothetical  3 inch  chord  faired  cable  of  NACA  0020  section  was  "towed" 
through  the  water  at  44  fps  (34.5  knots)  with  an  angle  of  attack  of  8°. 

Structurally,  the  towline  section  was  modelled  as  having  a rigid 

( ECORE  = 1x10^  psi)  CORE  extending  aft  from  the  leading  edge  for  3/4" 

( Y1 CORE  = 0.0,  Y2C0RE  = 0.75).  The  afterbody,  or  TAIL,  was  defined 

4 

to  be  a relatively  flexible  rubber-like  material  (ETAIL  = 1x10  psi). 

The  entire  fairing  section  was  enclosed  by  a very  flexible  SKIN 
(ESKIN  = 2x10^  psi),  0.05  inches  thick  (TSKIN  = 0.05).  The  finite 
element  mesh  used  to  analyze  the  structural  response  was  generated 
automatically  by  the  software  (MESHC  =11,  MESHT  = 3).  The  mesh 
consisted  of  48  2-0  elements  having  57  nodes;  the  outer  SKIN  connecting 
the  surface  nodes  was  modelled  by  24  1 -D  elements. 

For  the  purposes  of  the  hydrodynamic  computations,  the  fairing's 
NACA  0020  surface  profile  was  modelled  by  35  surface  points  distributed 
automatically  (HYDROSAP  controls  IPANEL  = 2,  and  N = - 3 ...  for  FOIL03) 
on  the  upper  and  lower  fairing  surfaces  by  the  software.  Viscous 
effects  were  included  in  the  computations  in  each  of  the  HYDROSAP 
"outer  iterations"  ( I TP,  SWT  = 1). 

The  effect  of  the  hydrodynamic  loading  on  the  flexible  symmetric 
fairing  was  to  induce  negative  camber.  Figure  5-1  summarizes  the  types 
of  deformations  suffered  by  the  fairing  cross-section.  Evidently,  the 
material  above  the  fairing  chord  was  compressed,  while  that  below  was 
in  tension.  The  equilibrium  trailing  edge  transverse  deflection  was 
approximately  2.3*  of  the  chord. 
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Figure  5-1.  TRANSVERSE  DEFLECTION  OF  CENTERLINE  FOR  A FLEXIBLE  CABLE  FAIRING 


The  negative  camber  acted  to  reduce  the  magnitude  of  the  lift 
and  (nose)  moment  coefficients  significantly.  In  Figure  5-2,  the 
hydrodynamic  coefficients  are  plotted  for  each  successive  pass  through 
the  HYDROSAP  system.  It  is  evident  from  the  graph  that  after  10 
iterations  the  hydrodynamic  loading  and  structural  response  are  in 
equilibrium.  The  fairing's  deformation  has  reduced  the  lift  coefficient 
by  24%,  from  = .8996  for  the  undeformed  fairing,  to  = .6794 
for  the  deformed  fairing.  The  nose-down  moment  coefficient,  taken 
about  the  leading  edge,  has  been  reduced  by  29%.  Figure  5-3  shows  the 
effect  that  the  modest  chordwise  deflection  of  the  fairing  has  on  the 
pressure  distribution.  The  pressure  coefficients  are  reduced  in  magni- 
tude by  about  10%  to  20%  even  in  the  rigid  leading  edge  region.  The 
trailing  edge  deflection  evidently  makes  its  presence  felt  globally  over 
the  fairing  surface. 


It  is  clear  that  the  chordwise  flexibility  of  a towline  fairing 
can  have  significant  effect  upon  the  fairing's  hydrodynamic  properties. 


FIGURE  5-2.  HYDRODYNAMIC  COEFFICIENTS  FOR  A FLEXIBLE  CABLE  FAIRING 
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FIGURE  5-3.  PRESSURE  DISTRIBUTIONS  OVER  UNDEFORMEO  ANO  DEFORMED  FAIRINGS 


6. 


References 


Stevens,  W.  A.,  S.  H.  Goradia,  andJ.A.  Braden;  "Mathematical  Model  for 
Two-Dimensional  Multi-Component  Airfoils  in  Viscous  Flow".  NASA 
CR-1843,  July  1971. 

Brune,  G.  W.  and  J.  W.  Manke:  "An  Improved  Version  of  the  NASA/ 
Lockheed  Mul ti -Element  Airfoil  Analvsis  Computer  Program".  NASA 
CR-145323,  March  1978. 

Morgan,  H.  L.,  Jr.:  "A  Computer  Program  for  the  Analysis  of  Multi- 
Element  Airfoils  in  Two-Dimensional  Subsonic,  Viscous  Flow". 

Aerodynamic  Analyses  Requiring  Advanced  Computers  Conference.  NASA/ 
Langley  Research  Center,  Hampton,  Virginia,  March  1975. 


(4)  Smetana,  F.  0.,  D.  C.  Summey,  N.  S.  Smith  and  R.  K.  Carden:  "Light 
Aircraft  Lift,  and  Drag,  and  Moment  Prediction  - A Review  and  Analysis". 
NASA  CR-2523,  May  1975. 

(5)  Oeller,  H.  J.:  "Die  inkompressible  Potential stromung  in  der  ebenen 
Gitterstufe  (Incompressible  Potential  Flow  in  Two-Dimensional  Cascades)". 
Wi ssenschaf tl iche  Gessel 1 schaf t fur  Luft-und  Raumfahrt  E.  V.,  pp. 

349-353,  Jahrbuch  1962  (In  German,  not  translated). 

(6)  Karamcheti , K. : Principles  of  Ideal-Fluid  Aerodynamics,  John  Wiley  & 

Sons,  Inc.,  New  York/London/Sidney,  1966. 

(7)  Cohen,  C.  B.  and  E.  Reshotko:  "The  Compressible  Laminar  Boundary  Layer 
with  Heat  Transfer  and  Arbitrary  Pressure  Gradient"  NASA  TND-1294,  1956. 

(8)  Schlichting,  H.:  Boundary  Layer  Theory.  Sixth  Ed.,  McGraw-Hill , New 
York,  1968. 


73 


(9)  Truckenbrodt,  E.:  "Ein  Quadraturuerfahren  zur  Berchnung  der  laminaren 
und  turbulenten  Reibungsschi ft  bei  ebener  und  rotationssymmetrischer 
Stromung".  Ing.  Arch.,  vol . 20,  pp  221-228,  1952  (NASA  TM-1379,  1955). 

(10)  Squire,  H.  B.  and  B.  A.  Young:  "The  Calculation  of  the  Profile  Drag 
of  Airfoils",  RAE  Report  R&M  No.  1838,  November  1937. 

(11)  K.  J.  Bathe,  H.  Ozdemir,  and  E.  L.  Wilson,  "Static  and  Dynamic 
Geometric  and  Material  Nonlinear  Analysis",  Report  No.  UC  SESM 
74-4,  Structural  Engineering  Laboratory,  Department  of  Civil 
Engineering,  University  of  California,  Berkeley,  California,  1974. 

(12)  K.  J.  Bathe,  E.  L.  Wilson,  R.  H.  Iding,  "NONSAP:  A Structural  Analysis 
Program  for  Static  and  Dynamic  Response  of  Nonlinear  Systems",  Report 
No.  UC  SESM  74-3,  Structural  Engineering  Laboratory,  Department  of 
Civil  Engineering,  University  of  California,  Berkeley,  California, 

1974. 

(13)  J.  T.  Oden,  Finite  Elements  of  Nonlinear  Continuua.  McGraw-Hill, 

1972. 

(14)  0.  C.  Zienkiewicz,  The  Finite  Element  Method  in  Engineering 
Science , McGraw-Hill,  London,  1971. 

(15)  Dvorak,  F.  and  B.  Maskew:  "The  Prediction  of  C using  a 

rri  Mnax 

Separated  Flow  Model".  Proceedings  of  the  33—  AHS  Meeting,  1977 
(Paper  No.  77.33.01.)  Also,  Journal  American  Helicopter  Society, 

April  1978. 


74 


APPENDIX  A 


HYDROSAP  Overlay  Structure 


HYDROSAP  Overlay  Structure 


In  order  to  retain  the  full  capabilities  of  the  HYDRO  and  NONSAP 
modules,  which  comprise  HYDROSAP,  and  yet  keep  the  entire  HYDROSAP  system' 
core  requirements  to  a minimum,  an  extensive  overlay  structure  has  been 
prepared.  Currently  the  system  requires  approximately  390  k bytes  of 
core.  A versior  without  overlay  would  require  over  600  k.  Maps  of 

the  overlay  structure  are  given  in  Figures  A-l,  A-2,  and  A-3. 


Figure  A-l.  HYDROSAP  Overlay  Tree 


Figure  A-2.  HYDRO  Overlay  Tree 


Figure  A-3.  NONSAP  Overlay  Tree 


APPENDIX  B 


HYDROSAP  Auxiliary  Subroutine  Descriptions 


Subroutine:  BEGIN® 


FUNCTION: 


BEGIN®  initializes  HYDROSAP  variables  and  sets  defau 
control  variables  brought  in  via  NAMELIST  / DRV/  (see 


CALLING  SEQUENCE: 

Call  BEGIN® 

COMMON  BLOCKS: 

/LOGIC0/ 
/REAL®/ 
/ARRAY®/ 

1/0  UNITS: 

None. 

PRINCIPAL  VARIABLES: 


/ INTGR0/ 
/ I NTGR1 / 
/1NTGR2/ 
/SURFCE/ 


See  subroutine  INPUT®. 


EXTERNAL  REFERENCES: 


None. 


CALLED  BY: 

DRIVER 
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t values  for  the 
subroutine  INPUT®). 


Subroutine  CONVRG 


FUNCTION: 

CONVRG  tests  the  size  of  the  current  displacement  increment,  relative  to 

the  total  displacement,  against  a user-assigned  tolerance  value. 

CONVRG  is  called  after  each  of  the  NONSAP  displacement  computations,  in 

each  HYDROSAP  outer  iteration. 

CALLING  SEQUENCE : 

Call  CONVRG  (DISP,  DISPI,  NEQ,  TOL,  1ST0P). 

COMMON  BLOCKS: 

/INTGR0/ 

/LOG1C0/ 

I/O  UNITS: 

Unit  IPAPER  = 6 for  printer  output 

= 22  for  dummy 

Unit  12  Terminal  summary  output. 

PRINCIPAL  VARIABLES: 

DISP  array  of  displacements  at  each  of  the  NONSAP  nodes 

(input) 

DISPI  array  of  displacement  increments  at  each  of  the  NONSAP 

nodes  (input) 

NEQ  number  of  equations  (active  degrees  of  freedom)  in 

NONSAP  (input) 

TOL  tolerance  level  assigned  for  the  convergence  test  (in- 

put) 
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I STOP  convergence  status  (output) 

= 0 not  converged;  proceed  with  next  HYDROSAP  outer 
iteration 

= 1 convergence  achieved;  normal  termination 
= 2 displacement  increments  all  zero;  error  termination. 

EXTERNAL  REFERENCES; 

*DSQRT  FORTRAN  Double-precision  square  root  function. 

CALLED  BY; 

DRIVER 


Subroutine  DEFALT 


FUNCTION: 

DEFALT  calls  the  appropriate  fairing  profile  library  routine  for  short 
form  input  ( SHRTI N= . TRUE . ) , based  upon  user  control  input  (variable  N). 

CALLING  SEQUENCE: 

Call  DEFALT  (X.Z.N). 

COMMON  BLOCKS: 

None. 

I/O  UNITS: 

None. 

PRINCIPAL  VARIABLES: 

X,Z  an  array  of  fairing  upper  surface  coordinates  expressed 

in  inches  (output). 

X(I)  = chordwise  location,  for  the  1th  station 
Z ( I ) = fairing  local  semi-thickness,  at  the  Ith  station 
(X(l ) , Z(1 ))  = leading  edge 
(X(N),  Z( N) ) = trailing  edge 
N = fairing  library  shape  routine  (input) 

= number  of  profile  points  (output). 

EXTERNAL  REFERENCES: 

FOIL01  Library  of  fairing  cross-section  (upper  profile) 

FO T L 02  arrays.  [Only  FOIL01  through  FOIL05  are  currently 

I implemented.] 
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FOIL05 


Program  DRIVER 

(alias  HYDROSAP) 

FUNCTION: 

DRIVER  is  the  main  HYDROSAP  executive  routine.  It  controls  the 
HYDRO  and  NONSAP  program  modules,  as  well  as  the  program  initialization 
input  file  preparation,  convergence  testing,  and  the  graphics  post- 
processor output. 


COMMON  BLOCKS: 


/INTGR0/ 

/LOGICA/ 

/ INTGR1 / 

/SURF CE/ 

/INTGR2/ 

/ARRAY#/ 

/REAL0/ 

/LOGIC#/ 

/ARRAY  1 / 

I/O  UNITS: 

Unit  5 

System  input  stream  ( SYS I N ) 

Unit  12 

Terminal  summary  output  (SYS0UT=X) 

Unit  IPAPER 

= 6 printer  output  (SYS0UT=A) 

= 22  dummy  output 

Unit  IREAD0 

NONSAP  Input  File  (=55) 

Unit  I READ! 

HYDRO  Input  File  (=66) 

PRINCIPAL  VARIABLES: 

SAPS  I M 

= .TRUE.,  NONSAP  simulation  mode  only 

HYDSIM 

=.TRUE.,  HYDRO  simulation  mode  only 

SHRTIN 

=.TRUE.,  short  form  input;  results  in  an  automatic  con 
struction  of  the  HYDRO  and  NONSAP  input  files 

TPRINT 

= .TRUE.,  terminal  print  summary  requested  ( SYSOUT  = X ) 

ECHO 

S.TRUE.,  "echo"  of  HYDRO  and  NONSAP  input  files  is 
produced 

POST 

=.TRUE.f  graphics  post-processor  output  is  produced 
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HYDROSAP  outer  iteration  counter  (value  written), 
maximum  allowable  value  for  ITER 
if  ITER.GE. ITRSWT,  then  the  HYDRO  computations  will 
include  viscous  effects;  otherwise  HYDRO  computations 
are  for  a potential  flow  solution. 

EXTERNAL  REFERENCES: 

BEGIN© 

1NPUT0 
HYOIN 
SAP  IN 
HYDRO 
NONSAP 
CONGRG 
OUTPST 

CALLED  BY: 

Not  applicable. 
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ITER 

LOOPMX 

ITRSWT 


Subroutine  FOILNN 


FUNCTION: 

FOILNN  is  a fairing  cross-section,  upper  surface,  library  routine  numbered 
"NN",  where  NN  can  range  from  01  to  99.  Files  denoted  "FOIL01"  through 
"FOIL05"  are  currently  implemented.  These  library  routines  will  sup- 
ply default  surface  coordinates  when  SHRTIN=.TRUE.  (see  DRIVER). 

CALLING  SEQUENCE: 

Call  FOILNN  (X,Z,N). 

COMMON  BLOCKS: 

None. 

I/O  UNITS: 

None. 

PRINCIPAL  VARIABLES: 

X,Z  arrays  giving  the  fairing's  upper  surface  coordinates, 

in  inches  (see  subroutine  DEFALT),  on  output. 

N number  of  fairing  upper  coordinate  pairs  on  the  output 

file. 

EXTERNAL  REFERENCES: 

None. 

CALLED  BY: 

DEFALT 
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FOIL  ROUTINES  CURRENTLY  IMPLEMENTED: 


FOIL01 

FO1L02 


FOIL03 

FOIL04 

FO1L05 


A typical  fairing  profile  (measured  from  the  actual  cable 
sample  given  to  BTS,  Inc.) 

Fairing  Model  I,  taken  from:  "Exp.  Det'n.  of  Hydrodynamic 
Loading  for  Ten  Cable  Fairing  Models",  by  R.  Folb, 

DTNSRDC  Rept.  4610,  November,  1975. 

NACA  0020  section  having  a 3 inch  chord 
NCS  0010  section,  having  a 3 inch  chord 
Reserved  for  program  expansion. 
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Subroutine  HYD1N 


FUNCTION: 

HYDIN  is  an  auxiliary  routine  which  builds  a complete  HYDRO  input  file, 
based  upon  values  for  the  control  variables  input  through  NAMELIST  /DRV/. 
HYDIN  is  activated  (only)  when  SHRTIN=.TRUE . 

CALLING  SEQUENCE: 

Call  HYDIN  ( I OUT , IPAPER,  X,Z,N,V,  ALPHA,  ECHO,  TITLE). 

COMMON  BLOCKS: 


None. 

I/O  UNITS: 

Unit  66  HYDRO  input  file 

Unit  6 Printer  output  (SYS0UT=A) 

Unit  22  Dummy. 

PRINCIPAL  VARIABLES: 

HYDRO  input  file  number 

output  file,  to  echo  the  HYDRO  input  file  being  built 
arrays  of  fairing  (upper  surface)  coordinates,  expressed 
in  inches  (see  DEFALT) . 

number  of  fairing  upper  surface  coordinates 
freestream  velocity  (ft/sec.) 
fairing  angle  of  attack  (degrees) 

=.TRUE.,  produce  an  "echo"  of  the  HYDRO  input  file 
80  character  (run)  title 


IOUT 

IPAPER 

X,Z 

N 

V 

ALPHA 

ECHO 

TITLE 


h 
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EXTERNAL  REFERENCES: 


square  root 

floating  point  conversion. 

CALLED  BY: 

DRIVER 


DEFALT 

DSQRT  FORTRAN  double  precision 

DFLOAT  FORTRAN  double  precision 


Subroutine  INPUT0 


FUNCTION: 

INPUT0  over-rides  the  default  values,  assigned  to  HYDROSAP  control 
variables  (see  BEGIN0),  with  control  values  read  in  through  NAMELIST 
/DRV/. 

CALLING  SEQUENCE: 

Call  INPUT0. 

COMMON  BLOCKS: 

/REAL0/ 

/ I NTGR0/ 

/ INTGR1 / 

1/0  UNITS: 

Unit  5 
Unit  6. 

PRINCIPAL  VARIABLES: 

Principal  variables  are  read  in,  via  namelist  format,  through  /DRV/. 

A full  description  of  the  control  variable  functions  and  defaults  is 
included  in  the  discussion  of  NAMELIST  /DRV/.  We  include  here  only  a 
list  of  the  relevant  variable  names: 


/INTGR2/ 
/LOG I C0/ 


NAMELIST  /DRV/ 

INTORD,  INTORX,  IPANEL,  IPAPER,  ISM'-,  ITRSWT,  IWRT0,  IWRT1,  LOOPMX, 
MESHC,  MESHT,  N,  NSTEPS, 

CORE,  DEBUG,  ECHO,  HYDSIM,  POST,  SAFSIM.  SHRTIN,  SKIN,  TAIL,  TPR1NT, 
ALPHA,  ECORE,  ESKIN,  ETAIL,  PCORE . PSK3N,  PTAIL,  PZERO,  RHO,  TOL, 
TSKIN,  Y 1 CORE , Y2C0RE. 
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AD-A071  202 


UNCLASSIFIED 

2*2  " 
AD 

AO  71 202 


BUSINESS  AND  TECHNOLOGICAL  SYSTEMS  INC  SEABROOK  MD  F/G 

PROGRAM  DEVELOPMENT  TO  STUDY  FAIRED  TOWLINES* (U) 

JAN  79  J B EADES • V MAJER  N00014-78-C-0410 


F/6  13/6 


BTS-IR-79-74 


ONR-CR298-003-1 


Subroutine  OUTPST 


FUNCTION: 

OUTPST  writes  the  pressure  distribution  and  shape  data  to  output  unit 
88;  there  it  can  be  ready  by  a graphics  post-processor  — such  as  pro- 
gram HSGRAF  --  in  a subsequent  job  or  job  step. 

CALLING  SEQUENCE: 

OUTPST. 

COMMON  BLOCKS: 

/INTGR0/ 

/REAL0/ 

/LOGIC0/ 

I/O  UNITS: 

Unit  88  output. 

PRINCIPAL  VARIABLES: 

REALM  character  array  for  (run)  title 


HYDROSAP  control  variables,  see  discussion  of 
NAMELIST  /DRV/ 


ZTITLE(20) 

HYDSIM 

SAPSIM 

SKIN 

CORE 

TAIL 


/SURFCE/ 
/ARRAY0/ 
/ARRAY! / 
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Subroutine  SAPIN 


FUNCTION: 

SAPIN  is  an  auxiliary  routine  which  builds  a complete  NONSAP  input 
file  based  upon  values  assigned  to  the  control  variables  (input) 
through  NAMELIST  /DRV/.  SAPIN  is  activated  if  SHRTIN= .TRUE. 

CALLING  SEQUENCE: 

Call  SAPIN  (IOUT,  IPAPER,  Y.Z.N,  ECHO,  TITLE,  NSTEPS,  SKIN,  CORE,  TAIL, 
MESHC,  MESHT,  ESKIN,  ECORE,  ETAIL , TSKIN,  PCORE,  PTAIL, 
PSKIN,  Y1CORE,  Y2CORE) . 

COMMON  BLOCKS: 

None. 

I/O  Units: 

Unit  55  NONSAP  input  unit 

Unit  6 Printer  output  (SYSOUT=A) 

Unit  22  Dummy 

PRINCIPAL  VARIABLES: 


IOUT 

IPAPER 

Y,Z 

N 

ECHO 

TITLE  (10) 


NONSAP  input  file  number 

output  unit  number  for  NONSAP  input  file  echo 
coordinate  arrays  describing  a fairing  upper  surface 
number  of  upper  surface  coordinate  pairs 
=.TRUE.,  echos  the  input  deck,  written  to  unit  IPAPER 
REAL*8,  80  character  (run)  title. 
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NSTEPS 

SKIN 

CORE 

TAIL 

MESHC 

MESHT 

ESK1N 

ECORE 

ETAIL 

TSKIN 

PCORE 

PTAIL 

PSKIN 

Y1CORE 

Y2COREj 


HYDROSAP  control  variables.  See  the  discussion  of 
NAMELIST  /DRV/  for  details. 


EXTERNAL  REFERENCES: 


None. 


CALLED  BY: 

DRIVER 
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Subroutine  TAINT 


FUNCTION: 

TAINT  is  a general-purpose  table  lookup  and  interpolation  algorithm. 
SOURCE  AND  AUTHOR: 

AMES  Research  Center  (anon.) 

CALLING  SEQUENCE: 

Call  TAINT  (XTAB,  FTAB,  X,  FX,  N,  K,  NER,  MON). 

COMMON  BLOCKS: 


None. 

I/O  UNITS: 

Unit  6 Printer  output  (SYSOUT=A) 

PRINCIPAL  VARIABLES: 

I I 

XTAB(l) 

FTAB(l) 

X 

FX 
N 
K 

NER 
MON 


1 

B-19 
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array  of  x-coordinates  (input) 
array  of  function  values  (input) 

x-coordinate  at  which  function  value  is  desired  (input) 

interpolated  function  value  (output) 

number  of  points  in  XTAB,  FTAB  (input) 

interpolation  order,  k = max  (N,9) 

error  condition  flag 

MON/O:  XTAB,  FTAB  unchanged  from  previous  call 
M0N=0:  New  XTAB,  FTAB  arrays. 


EXTERNAL  REFERENCES: 


None. 

CALLED  BY: 

HYDRO,  SAPIN,  LOADER 


I I 

I 

f ! 


1 
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Subroutine  LOADER 


FUNCTION: 

LOADER  controls  the  computation  of  a set  of  consistent  loads  which  are 
to  be  applied  to  the  foil  at  the  NONSAP  fairing  surface  nodes.  These 
loads  are  based  on  the  "computed  pressure  distribution"  from  HYDRO. 

CALLING  SEQUENCE: 

Call  LOADER  (NOD,  IDIRN,  NCUR,  FAC,  NLOAD,  P,  XCT,  ZCT,  NXUPTH,  NXLOTH) . 

COMMON  BLOCKS: 

/INTGR0/ 

/LOGIC®/ 

/SURFCE/ 

I/O  UNITS: 

None. 

PRINCIPAL  VARIABLES: 

N0D(1)  NOD(I)  = node  to  which  the  Ith  nodal  load  is  applied 

IDIRN(l)  IDIRN(I)  = direction  (2'y,  3'z)  for  the  1*^  nodal  load 

NCUR(l)  NCUR(I)  = time  function  (non-dimensional  load  profile) 

to  which  nodal  load  I is  referred 
FAC(l)  FAC ( I ) * value  of  nodal  load  I 

NLOAD  total  number  of  nodal  loads  being  applied 

P(l)  HYDRO-computed  pressure  distribution  array  (packed) 

XCT(l),  ZCT ( 1 ) HYDRO  surface  coordinate  arrays  (packed) 

NXUPTH(7,2)  1 upper  and  lower  surface  coordinate  index  storage  (needed 
NXLOTH (7,2)  ) to  unpack  XCT,  ZCT,  P) 

YUP(41 ) ,ZUP(41 ) NONSAP  upper  surface  coordinates 
YL0(41 ) ,ZL0(41 ) NONSAP  lower  surface  coordinates 
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EXTERNAL  REFERENCES 


U 


TAINT 

NDLLDS 

CALLED  BY: 

LOADS  (in  NONSAP  module). 


Subroutine  NDLLDS 


FUNCTION: 

NDLLDS  determines  the  consistent  set  of  nodal  loads  from  the  NONSAP 
surface  geometry  and  input  pressure  distribution. 

CALLING  SEQUENCE: 

Call  NDLLDS  (Y,Z,P,PY,PZ,N) 

COMMON  BLOCKS: 

None 

I/O  UNITS: 

None 

PRINCIPAL  VARIABLES: 

N number  of  NONSAP  surface  nodes 

Y,Z  surface  coordinate  arrays 

P pressure  distribution  array 

PY.PZ  consistent  nodal  load  arrays 

EXTERNAL  REFERENCES: 

None 

CALLED  BY: 


LOADER  (in  NONSAP  module). 


Subroutine  CPTOP 


FUNCTION: 

CPTOP  converts  the  non-dimensional  pressure  coefficients  into 
pressures  (for  loads  purposes). 

CALLING  SEQUENCE: 

Call  CPTOP  (CP,  P0,  RHO,  V,  P) 

COMMON  BLOCKS: 

None 

I/O  UNITS: 

None 

PRINCIPAL  VARIABLES: 


CP 

pressure  coefficient  (input) 

P0 

static  pressure  (psi)  (input) 

RHO 

fluid  density  (slugs/ft?)  (input) 

V 

fluid  freestream  velocity  (ft/sec) 

P 

stream  dynamic  pressure  (psi)  (output) 

EXTERNAL  REFERENCES: 

NONE 

CALLED  BY: 

HYDRO 
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Subroutine  WRT14 


FUNCTION: 

WRT14  is  a subprogram  used  to  write  the  HYDROSAP  checkpoint/restart 
tape,  FT14F001 . WRT14  concatenates  the  four  NONSAP  restart  tapes  on- 
to one  tape;  also,  it  adds  to  the  tape  the  HYDRO  checkpoint/restart 
data. 


CALLING  SEQUENCE: 

Call  WRT14  (TPRINT,  TIME,  IPAPER,  ITER,  XCT,  ZCT,  P,  NXUPTH,  NXLOTH) 


COMMON  BLOCKS: 


None 

I/O  UNITS: 


FT14F001 
FT 04 F 001 
FT08F001 
FT09F001 
FT13F001 


HYDROSAP  restart  tape 


NONSAP  restart  tapes 


PRINCIPAL  VARIABLES: 


TPRINT 

TIME 

IPAPER 

ITER 

XCT 

ZCT 

P 

NXUPTH 
NXLOTH  J 


terminal  print  control  flag  (logical) 
time,  at  the  checkpoint 
printer  output  unit  number 
HYDROSAP  iteration  count  at  checkpoint. 


HYDRO  data  for  checkpoint 
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EXTERNAL  REFERENCES: 


TPSCAN 

ERRSET 

CALLED  BY: 

NONSAP 


I 


1 

i 4 

k • 

SL 


routine  to  scan  unformatted  tapes  for  number  of  records 
and  record  length 

IBM  utility  sub-program  to  suppress  error  messages 
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Subroutine  TPSCAN 

FUNCTION: 

TPSCAN  is  a program  used  to  determine  the  number  of  4-byte  words  per 
record,  and  the  number  of  logical  records,  on  an  unformatted  I/O 
unit. 

CALLING  SEQUENCE: 

Call  TPSCAN  (IUNIT,  NREC , NWORDS,  IER) 

COMMON  BLOCKS: 

None 
I/O  UNITS: 

Tape  IUNIT,  determined  by  calling  sequence 

PRINCIPAL  VARIABLES: 

IUNIT 
NREC 

NWORDS (1000) 

IER 

EXTERNAL  REFERENCES: 

None 


CALLED  BY: 

WRT14 


unit  number  being  scanned 
number  of  records  on  IUNIT 

an  array  containing  the  number  of  4-byte  words  on  each 
record,  up  to  a maximum  of  1000  records 
error  return  code 


Subroutine  READ14 


FUNCTION: 

READ14  is  a program  to  unload  the  HYDROSAP  restart  tape  (FT14F00); 
READ14  is  called  when  RESTRT=.TRUE. 


CALLING  SEQUENCE: 

Call  READ14  (TPRINT,  TSTART,  ITER,  XCT,  ZCT,  P,  NXUPTH,  NXLOTH,  IER) 
COMMON  BLOCKS: 

None 

I/O  UNITS: 

FT14F001 
FT04F001 
FT08F001 
FT09F001 
FT13F001 


HYDROSAP  restart  tape 


NONSAP  restart  tapes 


PRINCIPAL  VARIABLES: 


IER 

TSTART 

TPRINT 

ITER 

XCT 

ZCT 

P 

NXUPTH 

NXLOTH 


error  return  code 
time  at  restart 

logical  terminal  print  control  flag 
HYDROSAP  iteration  count  at  restart 


HYDRO  data  for  restart 
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EXTERNAL  REFERENCES: 


None 
CALLED  BY: 
NONSPI 
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HYDROSAP  IBM  Job  Control  Language 
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HYDROSAP  IBM  Job  Control  Language 


The  job  control  statements  necessary  to  run  the  HYDROSAP  system  on  the 
NASA/Goddard  Space  Flight  Center's  IBM  360/91  machine,  are  discussed  in  this 
appendix.  These  control  statements  are,  of  course,  machine-dependent. 

A typical  job  consists  of  five  steps:  (1)  jobcard,  (2)  system  input 
stream  definition,  (3)  1 ink-overlay,  (4)  execution,  and  (5)  user  notification. 

Jobcard 


The  jobcard  contains  the  job  name,  accounting  information,  and  time 
estimates  for  the  run.  Changes  are  normally  made  only  to  the  job  name  and 
time  estimates: 


Example: 

//•YCSVMHjSfi  .inf:  c , T . hhhhhh . HO OHO!  . . 00*.  HOT  IFY*YC£YM 

jobname  account  number  time  estimate 

System  Input  Stream  Definition 


The  SYSIN  jobstep  defines  the  contents  of  the  system  input  stream, 
FORTRAN  unit  FT05F001 . The  input  stream  contains  card  images  of  the  input 

to  be  read  by  HYDROSAP.  The  contents  of  this  file  are  altered  for  each  new 
run,  as  described  in  section  2.4. 


Example: 


/ EXEC  SYSIN 
✓SYSIN  I' I'  ♦ 

NACA  0020  FAIRING  CROSS  SECTION:  SAMPLE  ANALYSIS 
I'.IiPV 

SHRT 1N®T  * I PR  I NT =5*  ITPSMT=1 . LOOPflXalS,  MESHC=1 1 * MESHT =S» 

N=— 3» NST EPS-20*  CORE =T * PEPNG=T • SK 1 N=T . TAI L=T  « TPP I NT  = T • 

T START  * 0 . It  0 • ALPHA** . D 0 • ECOPE  *!.  D 06  • E SK I N=£ . P OS  • ET  A I L = I . D 05  * 
PCORE® . 3H0  • PT  A I L* . 3H0 » PSK IN®.  3P0*  TOL  = I . D-3«  TSKIN*5.  P-02. 
'vFRSTP=22.  HO*  YICOPE*0.  DO*  Y2C0PE®.  75P0*  POST=F « CHKPNT  = T « 

I PANEL® 0. ISNTH=4 

I- END 


Link-Overlay 


The  LINK  jobstep  assembles  and  link  edits  the  various  program  segments, 
constituting  the  HYDROSAP  system,  into  a single  executable  load  module,  based 
on  the  linkage  editor  commands  which  define  the  overlay  structure.  Object 
modules  for  the  root  segment  and  two  overlay  segments  are  stored  in  load 
libraries  ZOFIMDRV,  ZOFIMHYD,  and  ZOFIMFEL,  respectively.  These  control 
statements  are  normally  not  altered. 


Example : 


-'Ll  Mr  EXEC  PGM*IEMl.PM?Me'l  I ST . MHP.  LET  • OVLY.  SI  2E=  t£40K , ?5\’y  ' . 
PEG  1 0N=c.'5i'iF  « C OND=  >’4  • LT  > 

'CYSPP1NT  til*  $YSQOT«H.  DCB»Bl KSIZE«3509.  CPHCE=  -TYL  « <3»  1 > - PLSE> 
--  ♦ : V SPP 1 NT  nit  SYSOUT  = P . I'C  B=BL  K$ I ZE  = 35  0 ? 

•VSYSLMOD  I'D  UN1T»?314.DISP=  -'NEW*  PftSSY  * DSNeSMivLODMOD  CGCFC>  » 

C F'hi'  E = >v.  l 44  • >: 4 i.i  • £ ft , i > • • . POUND'* 

•'SYSOTI  DD  ON  1 T «2*3 1 4 . <P«CE»  f 6 I 44 . (40*  £0*  O • • « POUND > 

•-  • CYCL1B  DD  DCN=CYC1  . DUMMY.  disp=$hp 

DUMMY . D 1 SP=SHP 
FOPTLIF..  D I SP=CHP 
FDPTLTB. D1 SP=SHP 
PL IL  I B«  D I S'P=S'HP 
FQPTSSP. DISP=SHP 
D S N= v YC  £ . LDhDL  ID.  HI  SF-CHP 
DCN=CYC3. I OHDL  1 B.  D1  SP=CHP 


/ / 

/ / 

-'■/L I B 
'SYSUIN 

INCLUDE 

INCLUDE 


DD 
DD 
DD 
D D 
DD 
DD 
DD 
DD 
DD 
DD  * 

L 1 B -TOF  IMDPVU 
LIB  •'ZOFIMHYD  ' 


d :n=:y: i 
DSN«SYSI 
DCN=SYS.-J 
dsn=sy:i 
ds n= - v : t 


INCLUDE 

INCEPT 

INCEPT 

OVEPLHY 

INCEPT 

incept 

OVEPLHY 
INCEPT 
INCEPT 
INCEPT 
1 NCEPT 
OVERLAY 
INSERT 
INCEPT 
OVERLAY 
INCEPT 
INCEPT 


LlBfZOFIMFEL'* 

MH I N . Th I NT . C ONVPG . OUT  PST 

PEHL  o.  JNTGPO. LOGIC  Cl.  INTGP1.  I NT  GP£ . CUF’F  C E * LOG  I CP  • HPPRY  0 
ONE 

HYDIN. ShPIN. DEFHLT. FOILOI • FOILdE*  FOILOC. FOIL 04. FOIL 05 

BEGINii.  INPUT 0 

ONE 

HYDPQ. INTEP 
CP  TOP 

PLL. RLPMCH. BLPYEP. BLK40. CFRLL • CF5* CNFLTP. GEOMT . HLM4 
LDhIu  PDTF  L 1.1  • POTUF . T E MP . T HE  T h J . TUP  DTP . V I SCDC « WORK 
Tl.lO 

H1PGM.  PEPDIT.  PPTGM.  GEOM.  PPTCF'D-  TPRNS.  PRNEL • CUPV 

SMOCP. GHPh. CUBSPL > 0 VP LHP 

Tl.lO 

PDTFL  . POTVL  • KI1TTH.  DCT EC'«  CHMBEP . SMOOTH.  THFVL  « OELLEP 
C I MC  OL  • COUPC'E  * COMPP 


(Continued) 
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OVERLAY 
INSERT 
INSERT 
OVERLAY 
1 NSERT 
OVERLAY 
INSERT 
INSERT 
INSERT 
INSERT 
INSERT 
INSERT 
INSERT 
OVERLAY 
I NSERT 
INSERT 
OVERLAY 
INSERT 
OVERLAY 
INSERT 
INSERT 
OVERLAY 
INSERT 
OVERLAY 
INSERT 
OVERLAY 
INSERT 
OVERLAY 
INSERT 


TWO 

BDLYR.  DEPV. LAMBL  - PROOT  . BLTRAN,  TIJRBL?  CONFBL  • DLIM 
CDNF5. C0NF7. CDNF® 

TMD 

LOADS  0 
ONE 

NONSAR. ERROR' EL CAL . ADDRES? ASSEM. LDAHEF . COLSOL , ECU IT 
NEWIiAV.  WRITE.  STRESS-  COLHT . ADDBAN.  ADDDI . ADDMA.  ADDCM 
MIJLT  • ELEMNT , ST  1 ME . TT I ME  < CAUCHY 
WRT14. TRSCAN 

EM.  MTMD3D.  VMISES,  DRRRA6*  DIMEL*  THREED*  TODIM*  D I SUER 
MAT  MOD . TARE I G . NSAR I . ELGFLG - DAMPNG • RRCON . CRACK . NORMS 
TAPES?  VAR.  JUNK,  LDA.  CONST . EL.  I'lM.  SOL*  X00001  * A 
THREE 

NONSP I . ORf'OEF  . I NRIJT . LOADS . NODMAS,  I N I TAL  - LDADER.  NDLLDS 

READ 14 

THREE 

TRUSS*  RUSS* MATBAR 
THREE 

TODMFE  • TDFE « I N I TWA.  MATPT£.  G'UADS . G'UADM.  DERI  Q,  FUNCT2 

MAXM I N • STSTL • RQ~ I NV . STSTN 

FOUR 

FLT£D3« IVTMOD" VTMOD 
FOUR 

ELT2D4. ICDMOD. CDMOD- DCRACK 
FOUR 

ELTcDC - IELRAL • ELRAL . MIDER 
FOUR 

ELTc'D? « I DRUCK  « DRUCK  . RRAGER 


OVERLAY  FOUR 

INSERT  ELTc’D®.  MOONEY 
OVERLAY  FOUR 

INSERT  ELTc'D®.  EL£I'l  0*  EL£I'  1 l.EL£Dl£ 

OVERLAY  THREE 

I NSERT  THPEDM . THDFE  • MATWRT  . I NTWA3  < 0UADS3  • QUADM3  - I'ER  I G'3 . FUNCT 

1 NSERT  STST3L • STST3N 
OVERLAY  FIVE 

INSERT  ELT3D3. ICM0D3. CM0D3D 
OVERLAY  FIVE 

INSERT  ELT3D5*  ELT3D6 
OVERLAY  THREE 

INSERT  FREG'S.  SECANT  . MLTRLY.  BANDET . RNOF'M.  WRTFRQ,  MAXAEX*  READ  ID 
ENTRY  MAIN 


Execution: 


The  GO  step  initiates  execution  of  the  complete  load  module  which  was 
produced  in  the  link  step.  The  GO  step  JCL  also  defines  the  FORTRAN  I/O 


units  which  will  be  accessed  by  the  program. 


Example: 

'"GO  EXEC  PGM=^.  LINK.  SYSLMDD*  PEG  I ON=380K«  CDND=  U-LT>' 

"GO.FTOIFOOI  DD  UNI T=£31 4. SPPCE=  CTRK - 1 6 0 • 6 OV  :>  . 

/v  DCB=  (RECFM=VSB.  LPECL=  1 284 « BLKS  I ZE=257£>  - V0L=SER=M£SCR1 
. • - GO . FT  02F  0 0 1 DD  LIN  I T=£3 1 4 , SPACE  = <TRK  • <6  0 • 6 O'1'  > - 
/v  DCB=  <PECFM=VSB » LPECL  = 1 884 , BLKS I ZE=£57£>  . V0L=SER=M£SCR2 
•'•'GO.  FT03F001  DD  UNIT =2314.  SPPCE=  CTRK • (60*  GOV > . 

••'V  DCB=  <RECFM=VSB*  LPECL  = 1284.  BLKS I ZE=257£>' , V0L=SER=M£SCR3 
•'•'GO.  FT04F001  DD  UNIT  =2314 « SPACE®  f TPK « (60«60?>'» 

/ v DCB=  fRECFM=VSB • LPECL=  1 284 . BLKS I ZE=£572>  , V0L=SER=M2SCR4 
- "'GO . FT  05F  0 0 1 DD  DDN8ME=D8T PS 

••■  ••-GO . FT  08F  0 0 1 DD  SYS0UT=8 . DCB=  •'PECFM=VB8  • LPECL  = 1 37 , BLKS  I ZE=7£t5> 

-•'  • GO . F T 07 F 0 0 1 DD  DUMMY 

• 'GO.  FT08F001  DD  UN  I T *23 14*  SPACE*  * TRK » >'GO«GO>  • « 

DCB= <RECFM=VSB * LPECL= I £84 , BL KS I ZE=£57£> . VOL=SER=M£SCR£ 

"GO. FT 09F 001  DD  UNIT=831 4. SP8CF= (TRK.  (80«G0>  '. 

■ ■•■•  DCB=  ■:PECFM=VSB . L RECL  = 1 284 . BLKS  I ZE  =£572  > . OL  = S ER  = M2SCR3 

• GO. FT  I OF  001  DD  ONI T=231 4. S’PPCE  = fTPK * GO1  1 • 

••'■•••  DCB=  fRECFM=VSB<  LPECL  = 1£84.  BLFSIZE=2S.7£  * . v,0L»SER*N2SCR4 
'GO.  FT!  I Foul  DD  UNI  T=2314»  S:F'PCE=  ■' TRK.  •’►.i.i*  GO  1 '■  * 

DCB=  <rRECFM=VSB.  LRFCL  = 12&4.  BLKS  I ZE=2S~2>  . VOL  = TER=M2SCR1 
•GD.  FT  1 £F 001  DD  SYSDUT=X.  DCB=  fRE C FM=FB . L REC  L =80 • BLKS  IZE*8CO 
•■'GO . FT  1 3:F 0 0 1 DD  UNIT— £2:14«  S-PPCE=  *’TRK « '’GO.GO*  1 « 

•■'••'  DCB=  (RECFM=VSB.  LRECL=1284.  BLKS.1  ZE=257£  • . VDL*SER=M£SCR1 
♦GO.  FT  I 4F  in'il  DD  IJNIT=c‘3l4i  SPACE = '’TRC < ■' G 0 ■ G 0 1 > • 

♦ DCB=  (RECFM=VSB. LPECL =1284 . BLKS  I 2E*25?£  • . VOL  =SER=M£SCR£ 

•'GO.  FT14FOOJ  DD  DSN=YC2VM . US  1 4 . D8T8 • UN  1 T=£3 1 4 • SP8CE=  -'TRK.  <£0«20  • 
H ISP*  <NEW.PASS>  « DCB*  <RECFM«VSB«  LPECL*  1 £84.  BLKSIZE=£57£'-  . 

•v  V0L=SER=M£SCR2 

• GO . FT22F  0 0 1 DD  DUMMY 

■■'GO . FT55F  001  DD  UN  1 T =£'■'•  14.  SPACE  = >'  T PK  • f£0*20V.. 

•■  DCB*  ,’R£CF  M=FB . L RECL  =8  0 . ELKS  I ZE  =32  0 0 > • VOL  =S  ER=M£SCR£ 

-GO. FTGGF 001  DD  UNI T=231 4«  SPACE  = f TRK « f£0*  £0V>  . 

•"■  D C B = «•’  R E CF  M = F B « L P E C L = 8 0 « B L K S 1 ZE = 3 £ 0 0 • - V 0 L = S E R = M2  SC  R £ 

• ♦GO.  FT88F001  DD  UNIT  =£31 4*  SPACE*  »’TRK . .'80-  £0  ‘ > « 

•-•♦DCB*  (PECFM*FB«  LPECL  =80 • BLKS  I ZE  = 3£00  • • VOL  = SER=M£SCR3 

'GO.  FT88F001  DD  DSN=YC2VM.  HS88.  DAT 8.  UNI  T = £31 4.  SPACE*  ,rTPK « >'£0.  20 V 
DI  SF  = f NEW.  PAS S • • DC B*  ''RECFM=FB«  LRECL=80*  BLKS  IZE=3£00  - • VOL  = SER=M 
GO.  DPT 85  DD  DSNs&jj-DPT 85*  DISP*  • OLD-  DELETE'* 
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User  Notification 


The  NOTIFY  step  is  used  to  send  a notification  to  the  user's  remote 
terminal  that  the  program  has  completed  execution.  The  system  completion 
code  for  each  job  step  is  displayed. 

Example: 


EXEC  NOT  I FYTS « I Ii=YC2VM * MOI‘E=*LL  • " HYPROSfiF'  S-ftMPi-E 


« I ft 
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